arXiv: 1503.07530vl [cond-mat.mes-hall] 25 Mar 2015 


Proximity Induced Vortices and Long-Range Triplet Supercurrents in Ferromagnetic Josephson 

Junctions and Spin Valves 


Mohammad Alidoust^’^ * and Klaus Halterman^’ 

' Department of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, Switzerland 
^Department of Physics, Faculty of Sciences, University of Isfahan, Hezar Jerib Avenue, Isfahan 8I746-7344I, Iran 
^Michelson Lab, Physics Division, Naval Air Warfare Center, China Lake, California 93555, USA 

(Dated: March 27, 2015) 

Using a spin-parameterized quasiclassical Keldysh-Usadel technique, we theoretically study supercurrent 
transport in several types of diffusive ferromagnetic(F)/superconducting(5) configurations with differing mag¬ 
netization textures. We separate out the even- and odd-frequency components of the supercurrent within the 
low proximity limit and identify the relative contributions from the singlet and triplet channels. We first con¬ 
sider inhomogeneous one-dimensional Josephson structures consisting of a uniform bilayer magnetic S/FIF/S 
structure and a trilayer S/F/F/F/S configuration, in which case the outer F layers can have either a uniform 
or conical texture relative to the central uniform F layer. Our results demonstrate that for supercurrents flowing 
perpendicular to the FIF interfaces, incorporating a conical texture yields the most effective way to observe 
the signatures of the long-ranged spin-triplet supercurrents. We also consider three different types of finite¬ 
sized two-dimensional magnetic structures subjected to an applied magnetic field normal to the junction plane: 
a S/F/S junction with uniform magnetization texture, and two S/F/F/S configurations with differing F/F 
bilayer arrangements. In one case, the F/F interface is parallel with the S/F junction interfaces while in the 
other case, the F/F junction is oriented perpendicular to the S/F interfaces. We then discuss the proximity 
vortices and corresponding spatial maps of currents inside the junctions. For the uniform S/F/S junction, we 
analytically calculate the magnetic field induced supercurrent and pair potential in both the narrow and wide 
junction regimes, thus providing insight into the variations in the Fraunhofer diffraction patterns and proximity 
vortices when transitioning from a wide junction to a narrow one. Our extensive computations demonstrate that 
the induced long-range spin-triplet supercurrents can deeply penetrate uniform F/F bilayers when spin-singlet 
supercurrents flow parallel to the F/F interfaces. This is in stark contrast to configurations where a spin-singlet 
supercurrent flows perpendicular to the F/F interfaces. We pinpoint the origin of the induced triplet and singlet 
correlations through spatial profiles of the decomposed total supercurrents. We find that the penetration of the 
long-range spin-triplet supercurrents associated with supercurrents flowing parallel to the F/F interfaces, are 
more pronounced when the thickness of the F strips are unequal. Lastly, if one of the S terminals is replaced 
with a finite-sized normal metal, we demonstrate that the corresponding experimentally accessible S/F/F/N 
spin valve presents an effective platform in which the predicted long-range effects can be effectively generated 
and probed. 

PACS numbers: 74.50.+r, 74.25.Ha, 74.78.Na, 74.50.-l-r, 74.45.+C, 74.78.FK, 72.80.Vp, 68.65.Pq, 81.05.ue 


I. INTRODUCTION 

The interaction between the different order parameters in 
proximity coupled nanostructures comprised of ferromagnets 
(F) and superconductors {S) has attracted considerable atten¬ 
tion from numerous scientific disciplines in both the theoret¬ 
ical and experimental communities.*'® The interplay between 
ferromagnetism and superconductivity at low temperatures 
has constituted a unique arena for researchers in condensed 
matter studying superconducting spintronics in the clean, dif¬ 
fusive and non-equilibrium regimes.***'^* Interest in supercon¬ 
ducting electronics involving S/F hybrids has substantially 
increased during the past decade due to considerable advances 
in nanofabrication techniques. This consequently has led to 
more possibilities for S/F heterostructures playing a prac¬ 
tical role in nanoscale systems including, quantum comput¬ 
ers and ultra-sensitive detectors.Several interest¬ 
ing and important effects have been found and studied both 
theoretically and experimentally, such as O-tt transitions^®'^^, 
and the existence of triplet correlations"*’^’^^'^*. 

When a single quantization axis can be defined through¬ 
out the system, such as in a simple S/F bilayer with a 


uniformly magnetized ferromagnet, the Cooper pair wave- 
function is composed of singlet and opposite-spin triplet 
components."*’^^’^"* These components have zero spin projec¬ 
tion along the quantization axis, which is the same direc¬ 
tion as the magnetization. These two types of supercon¬ 
ducting correlations oscillate and strongly decay inside the 
F layer over length scales determined by fp. In the diffu¬ 
sive regime studied here, fp — ^jD/h, where D and h rep¬ 
resent the diffusion constant and exchange field magnitudes, 
respectively. In the ballistic regime, /,p = hvp/2h, where 
vp is the fermi velocity. Due to these relatively small length 
scales, the zero-spin triplet correlations are often referred to 
as short-ranged."*’^’^^’^® However, if the magnetization of the F 
layer possesses an inhomogeneous pattern, equal-spin triplet 
correlations can be generated.These types of corre¬ 
lations have non-zero spin projection along the quantization 
axis. The equal-spin correlations penetrate into a uniform dif¬ 
fusive F media over a length scale that is the same as singlets 
in a normal metal.For instance, it has been theoreti¬ 
cally shown that in the diffusive regime, a particular trilayer 
S/F/F/F/S Josephson junctions with non-collinear magne¬ 
tizations may host triplet supercurrents that are manifested 
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in a slowly decaying critical current as a function of junc¬ 
tion thickness^^. It has also been demonstrated that to reveal 
the long-ranged nature of proximity triplet supercurrents in 
the diffusive regime, a simple uniformly magnetized S/F/F/S 
junction may not possess the requisite magnetic inhomogene¬ 
ity, and consequently a counterpart layered SIFIFIFIS junc¬ 
tion is necessary"*^. In contrast to the diffusive regime, it has 
recently been shown in the ballistic regime that it is possi¬ 
ble to generate long-range odd-frequency triplet correlations 
in SIF/FIS Josephson junctions containing two uniform F 
layers with misaligned magnetization orientations and differ¬ 
ing thicknesses [dpi 7^ The signature of these 

long-ranged triplet correlations are theoretically predicted to 
be revealed in the second harmonic term of the Josephson 
current'^^’"'^’'^^. The signatures of the equal-spin triplet corre¬ 
lations have been observed in experiments as 
When the magnet is fully spin polarized, as in half-metallic 
systems, these type of triplet correlations can be produced 
when there are spin-active interfaces present.^*’^^’^'’^^ A sig¬ 
nificant thrust of these works is the formulation of simple 
and optimal conditions to detect the odd-frequency pairings 
in S/F heterostructures. To this end, spin-valve S/F/F het¬ 
erostructures have recently attracted interest from both the 
theoretical and experimental communities^®’'*^’^^’^^^®®. The ad¬ 
vantages of such spin valves are their less complicated experi¬ 
mental implementation and greater control of their magnetiza¬ 
tion state compared to layered magnetic Josephson junctions. 

In this paper, we make use of a spin parametrization scheme 
for the Green’s function, the Usadel equation, and associ¬ 
ated boundary conditions. This method provides a suitable 
framework for separating the supercurrent into spin singlet, 
opposite-spin triplet, and equal-spin triplet components, us¬ 
ing the spin parametrization technique outlined for a generic 
three dimensional system. Our model allows for investi¬ 
gations into a broad range of realistic finite-size ferromag- 
net/superconductor hybrids with arbitrary magnetization pat¬ 
terns subject to an external magnetic field^’^^’^^. The spin de¬ 
composed supercurrent accurately pinpoints the contribution 
from different superconducting pairings, their influence upon 
the total supercurrent, and their spatial variations within the 
magnetic regions.^ 

We first consider three types of one-dimensional Joseph¬ 
son junctions and study the critical supercurrent spin decom¬ 
posed components for differing F layer thicknesses. Our 
results demonstrate that in the low proximity regime of the 
diffusive limit, the most effective way to observe signatures 
of the long-ranged spin-triplet supercurrents (where the su¬ 
percurrent flows perpendicular to F/F interfaces) involves 
the use of inhomogeneous magnetic structures comprised of 
combinations of rotating exchange interactions (e.g. conical 
texture in Holmium [Ho]) and uniform ferromagnets. We 
find that the supercurrent spin decomposed component corre¬ 
sponding to the rotating component of magnetization is long- 
ranged in such a situation and dominates the behavior of to¬ 
tal supercurrent. Trilayer SIFIFIFIS structures with uni¬ 
form ferromagnets"^^’®^ are shown to weakly display long- 
range spin-triplet signatures in this low proximity limit. 

Next, we consider three different types of finite-sized fwo- 


dimensional magnetic Josephson junctions subject to an ap¬ 
plied magnetic field. ^ Our general analytical and numerical 
framework permits the study of magnetization textures with 
highly intricate patterns.Our methodology also allows 
for rather general geometric parameters, including arbitrary 
ratios of the side lengths describing the ferromagnet strips. 
We first consider a SI FIS Josephson junction with a uniform 
magnetization texture, thus extending the results of a normal 
SINIS Josephson junction.®^ In doing this, we employ simpli¬ 
fying approximations that permit explicit analytical solutions 
to the anomalous Green’s function. This consequently leads 
to tractable and transparent analytical expressions for the spa¬ 
tial dependence to the current density and pair potential. In 
particular, we implement the so-called wide and narrow junc¬ 
tion limits, which results in considerable simplifications to the 
Usadel equations. In the wide-junction limit, the Fraunhofer 
diffraction pattern appears with $0 (the magnetic flux quan¬ 
tum) periodicity in the critical supercurrent as a function of 
external magnetic flux, whereas a narrow-junction transitions 
from an oscillating Fraunhofer pattern to a monotonically de¬ 
caying one, similar to its normal metal counter part®*. As¬ 
sociated with these signatures of the supercurrent is the ap¬ 
pearance of arrays of proximity vortices®*’^'^^"^ which provide 
useful information regarding the Fraunhofer response of the 
supercurrent to an external magnetic held. By initially consid¬ 
ering uniformly magnetized structures with a single F layer, 
the nature of the proximity vortices and current flow mappings 
in more complicated magnetically inhomogeneous junctions 
discussed below are better understood in addition to the O-tt 
transition influences on the critical supercurrent responses^”. 

To explore the possibility of induced long-range triplet ef¬ 
fects, additional magnetic inhomogeneity is introduced by the 
addition of another ferromagnet layer, thus establishing dou¬ 
ble magnet SIFIFIS Josephson junctions. These types of 
structures comprise the main focus of the paper. Two types of 
SIFIFIS configurations are considered: In one case, the FIF 
interface is parallel to the interfaces of the S leads, while in the 
other case, the FIF junction is oriented perpendicular to them. 
In either scenario, when an external magnetic held is present, 
it is applied normal to the junction plane. Our findings demon¬ 
strate that a diffusive SIFIFIS Josephson junction in the low 
proximity limit can generate long-ranged triplet supercurrents 
depending on the direction of charge supercurrent with respect 
to the FIF interface orientation. In particular, if charge su¬ 
percurrent flows parallel with the FIF interface, spin-triplet 
components generated in one ferromagnetic wire deeply pen¬ 
etrate the adjacent ferromagnet with relative orthogonal mag¬ 
netizations. For these types of structures, we And that the 
long-ranged effect manifests itself when the thickness of the 
ferromagnetic strips are unequal.^ With the goal of demon¬ 
strating the generality of the introduced scenario above, iso¬ 
lating the predicted equal-spin triplet component to the su¬ 
percurrents flowing parallel to FIF interfaces, and motivated 
by recent experiments involving S/F/F spin-valves^®’®'’®®^®, 
we turn our attention to S/F/F/N spin-valves subject to an 
external magnetic held {N denotes a normal metal layer). Our 
results show that indeed for certain geometric and material pa¬ 
rameters, diffusive S/F/F/N spin valves can isolate purely 
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equal-spin odd-triplet correlations arising from the Meissner 
response, following the parallel transport scenario above, even 
in the low proximity limit. The supercurrent moving parallel 
to FIF contact, in this case is long-ranged, extends consid¬ 
erably into the N layer, and can be experimentally probed 
through direct local measurements of the current inside the 
relatively thick normal layer. We find that an equal-spin triplet 
supercurrent appears when the thickness of the two F layers 
are unequal: dpi ^ dF 2 , and vanishes when dpi = dp 2 , 
consistent with the behavior of the ferromagnet Josephson 
junctions mentioned above. Therefore, our extensive study 
demonstrates the generality of our proposed scenario to ef¬ 
fectively generate long-ranged supercurrents independent of 
geometry implemented. 


The paper is organized as follows: We present a succinct 
review of the theoretical framework, spin-parametrization 
scheme, and parameters employed in Sec. 11. In Sec. Ill, the 
one-dimensional spin-parameterized Green’s function is dis¬ 
cussed, and in Sec. III. 1 we present the approach taken to 
evaluate the corresponding decomposed supercurrents. In 
Sec. III.l, we discuss the critical supercurrent, O-tt transi¬ 
tions and equal-spin triplet components of the supercurrent for 
one-dimensional S/F/FIS, S/F/F/F/S, and S/Ho/F/Flo/S 
structures. Next, in Sec. IV we expand our investigations into 
two-dimensional hybrid junctions. In Sec. IV. 1, the pertinent 
technical points and parameters used to study the proposed 
heterostructures theoretically are presented, and which are 
chosen to be aligned with realistic experimental conditions. 
In Sec. IV.2, we consider i) the wide-junction limit, and ii) 
the narrow-junction limit of a uniform S/F/S junction. Cor¬ 
responding analytical expressions are given for the pair poten¬ 
tial and supercurrent response when the system is subject to an 
applied magnetic field. In part iii) we compliment our analyt¬ 
ical expressions with a full numerical treatment that does not 
resort to the previous simplifying assumptions. In Sec. IV.3, 
we present the spin-parameterized Usadel equation, supple¬ 
mentary boundary conditions, and separate out the contribu¬ 
tions from the odd and even frequency components of the net 
supercurrent describing these two-dimensional systems. In 
Sec. IV.4, we study one of the main structures, a magnetic 
S/F/FIS junction, where the double layer F/F interfaces are 
aligned with the interfaces of the S terminals. In Sec. IV.5, 
the remaining structure is discussed, where the F/F interfaces 
are orthogonal to the interfaces of the S terminals. For both 
configurations, we study the pair potential, charge supercur¬ 
rent, and its odd or even frequency decomposition. We show 
also how an external magnetic flux can induce vortex phenom¬ 
ena and modifications to the singlet and triplet correlations 
responsible for supercurrent transport. We also study the in¬ 
fluence of ferromagnetic strip thicknesses on the long-range 
spin-triplet contributions to the charge supercurrent. Finally, 
in Sec. IV.6, we study the long-range spin-triplet supercur¬ 
rents in S/F/F/N valves. These results are then compared 
with those obtained for their SIFIF/S counterparts. The ex¬ 
perimental implications of our findings for these structures are 
also discussed. Finally, we summarize our findings in Sec. V 
and give concluding remarks. 


II. GENERAL APPROACH AND EORMALISM 

Here we first outline the theoretical approach for generic 
three-dimensional systems. The corresponding reduced one¬ 
dimensional and two-dimensional cases are presented in the 
subsequent sections. 


II. 1. Theoretical methods 


The coupling between an s-wave superconductor and a 
ferromagnet leads to proximity-induced triplet correlations 
in addition to the usual singlet pairings.The correspond¬ 
ing coherent superconducting quasiparticles inside a diffusive 
medium can be described by the Usadel equations,which 
are a set of coupled complex three-dimensional partial dif¬ 
ferential equations. A general three-dimensional quasiclassi- 
cal model for such diffusive feiromagnet/superconductor het¬ 
erostructures subject to an external magnetic held is given by 
the following Usadel equation;'*’^^’^® 

D[d,G{x,y, z)[d,G{x,y, z)]] + i[ep 3 + 

dmg[h{x, y, z) ■ a, {h{x, y, z) ■ a)'^],G(x, y, z)] = 0, (1) 

in which and a are 4x4 and 2x2 Pauli matrices, 
respectively, and we denote the diffusive constant of the 
medium by D. Here, the exchange held of a ferromagnetic 
region, h{x,y,z) = {h^{x,y, z),hy{x,y, z),h^{x,y, z)), 
can take arbitrary directions in configuration space. We 
have defined the 4x4 version of partial derivative, d, by 
d = VI — ieA{x,y, z)p 3 in which A{x,y,z) stands 
for a vector potential producing the applied external mag¬ 
netic held H{x,y,z) and [d,G[x,y,z)\ = yG{x,y, z) — 
ie[A{x,y, z)p 3 ,G{x,y, z)]. We have denoted the quasipar¬ 
ticles’ energy by e which is measured from the fermi surface 

£p. 

In the low proximity limit, the normal and anomalous com¬ 
ponents of the Green’s function can be approximately written 
as, lS^°{x,y,z) ~ 1 and F{x,y,z) ^ 1, respectively. In 
this regime therefore, the advanced component of the Green’s 
function can be directly expressed as: 


G^{x,y,z,e) 


-1 -F{x,y,z,-£) 

F*{x,y,z,£) 1 


( 2 ) 


where the underline notation reflects 2x2 matrices. Thus, the 
advanced component, G^{x,y, z), of total Green’s function 
G can be written as: 


G^{x,y,z) 


/ -1 0 -/n(-£) -/-(-e)\ 

0 -1 -/+(-£) -/u(-e) 

/f(e) flie) 1 0 

\.ni^) 0 1 


In general, when a system is in a nonequilibrium state, the 
Usadel equation must be supplemented by the appropriate 
distribution functions.In this paper, however, we assume 




4 


equilibrium conditions for our systems under consideration, 
and hence the three blocks comprising the total Green’s 
function are related to each other in the following way: 

G^ix,y,z\ = -{T 3 G^{x,y,z)T 3 )K and G^{x,y,z) = 
ta,iih{j3e){G^{x,y, z) — G^ {x, y, z) ), where ra is the Pauli 
matrix, and (3 = kBT/2. 

The resulting nonlinear complex partial differential equa¬ 
tions should be supplemented by appropriate boundary 
conditions to properly capture the electronic and trans¬ 
port characteristics of SIF hybrid structures. We employ 
the Kupriyanov-Lukichev boundary conditions at the SIF 
interfaces^* and control the induced proximity correlations us¬ 
ing the parameter as the barrier resistance: 

C(G(a;, y, z)dG(x, y, z)) ■ h = [G'bcs(6'), G{x, y, z)] + 
i(Gs/Gt)[ diag(73, ^), G{x, y, z)], (4) 

where n is a unit vector denoting the perpendicular direction 
to an interface. The parameters Gg and Gt introduce spin- 
activity at the SIF interfaces.The solution for a bulk even- 
frequency s-wave superconductor G|cs reads, 

rR (a\-( lcoshr?(e) ir 2 e*® sinh■(?(£) \ 

Bcsl ) sinhi9(£) —lcosh')9(£) J ’ ^ 

where, 

Die) = arctanh(J-^). 

e 

Here we have represented the macroscopic phase of the bulk 
superconductor by 6. To have more compact expressions, we 
define the following piecewise functions: 


Wp is smaller than the Josephson penetration length 
Therefore, screening of the magnetic field by Josephson cur¬ 
rents can be safely ignored.In general, the external magnetic 
field strongly influences the macroscopic phases of the super¬ 
conducting leads. To avoid such effects in our considered sys¬ 
tems, we assume that the external magnetic field passes only 
across the nonsuperconducting sandwiched strips. These as¬ 
sumptions lead to results that are in very good agreement with 
those found in experiments.This will be discussed 
in more detail in Sec. IV. 

One of the most important quantities in the context of 
quantum transport through Josephson junction systems is 
the charge supercurrent, which provides valuable information 
about the superconducting properties of the system and the 
associated favorable experimental conditions under which to 
detect them. Under equilibrium conditions, the vector current 
density can be expressed by the Keldysh block as follows: 

/ +00 

deTrlps {G{x, y, z) [d, G{x, y, z)])^ |, 

-OO ^ 

(6) 

where Jg = NqcD/A, Nq is the number of states at the Fermi 
surface, and e is the electron charge. The vector current den¬ 
sity, J{x,y,z), provides a local spatial map and measure of 
the charge supercurrent flow through the system. To obtain 
the total Josephson charge current flowing along a particular 
direction inside the junction, it is necessary to perform an ad¬ 
ditional integration of Eq. (6) over the direction perpendicular 
to the transport direction. For instance, the total charge current 
flowing along the x direction can be obtained from, 

Ix=Io dyj{x,y,z) ■ X, (7) 

JWf 


sjs) = sinh ■!?(£) = 




i 

c(£) = coshi?{£) = 

I e ' 


\/£^ — 


0(£^ - A") - 


le 


\/A^ 


=0(A^-e^), 


where 0(a:) denotes the usual step function. 

In the situations where an external magnetic held is applied, 
it is directed along the z-axis. We also use the Coulomb gauge 
V • A{x, y,z) = 0 throughout our calculations for the vec¬ 
tor potential. In a magnetic junction the vector potential is 
composed of two parts: a) a part due to the magnetic held 
associated with the exchange interaction in the ferromagnetic 
layer h{x, y, z), and b) a part due to the external magnetic 
held H{x,y,z). The contribution due to the exchange in¬ 
teraction measured in experiments reveals itself as a shift in 
the observed magnetic interference patterns^However, 
as this has found good agreement with experiments, one can 
safely neglect the part of A{x, y, z) arising from the exchange 
interaction.Thus, we assume that the external magnetic 
flux contribution dominates, and the vector potential can be 
determined entirely by the external magnetic held H{x^y,z). 
In this paper, we consider the regime where the junction width 


in which Wp is the junction width (see e.g.. Fig. 5). Likewise, 
the charge supercurrent flow in the y direction can be obtained 
via integration of Jy{x, y, z) over the x coordinate. Another 
physically relevant quantity which gives additional insight 
into the local behavior of the singlet correlations throughout 
the Josephson structure is the spatial maps of pair potential, 
Upmr{x, y, z). This pair correlation function is defined using 
the Keldysh block of the total Green’s function, G{x, y, z)^®; 



( 8 ) 


where we normalize the pair potential as, Upaaix, y, z) = 
Upaii{x,y, z)/Uo, where Uq = —NqX/S, and A is a constant 
inside the superconducting regions. We note that although the 
pair potential must vanish outside of the intrinsically super¬ 
conducting regions, the pair amplitude, Up-^nix, y, z), is gener¬ 
ally nonzero in the ferromagnetic regions due to the proximity 
effect. 


II.2. Spin-parametrization and parameters 

Due to the possible appearance of triplet pairings in such 
hybrid structures^, we assume a fixed quantization axis and 












5 



FIG. 1. (Color online) Schematic of three experimental sets 
up. (a) simple S/F/F/S Josephson junction with ferromagnetic 
layers’ widths dpi, dF2 and magnetization orientations hi = 
/io(0, sin/3i, cos/3i), /12 = ho(0,0,1), respectively, (b) trilayer 
S/F/F/F/S junction where the magnetization of middle F layer, 
with width dF 2 , is fixed in the 2 direction. The other layers 
have dpi, dps widths with magnetization orientations = 

/io(0, sin/3i,3, cos/?i, 3 ), respectively, (c) trilayer of ferromagnets 
where the outer F layers have Holmium-like magnetization patterns 
namely, hno ~ /io(cos a, sin a sinyoj/a, sin a cos 7 a:/a). The 
length of middle F layer is labeled hy dp while the two outer 
Holmium-like F layers are assumed to have equal lengths dpo- The 
X axis is normal to the junction interfaces and the systems have infi¬ 
nite sizes in the y and z directions. 


employ a spin-parametrization scheme. In this scheme, the 
Green’s function is decomposed into the even and odd fre¬ 
quency components, taking the spin-quantization axis to be 
oriented along z direction. Since we consider the low prox¬ 
imity limit, the anomalous component of the Green’s function 
takes the following form in terms of even- (S) and odd- (T) 
frequency parts: 


y, z,e) =i S(x, y,z,e)+T- T(a;, y, z, e) 


(9) 


where r = {Tx,Ty,Tz) is a vector comprised of Pauli ma¬ 
trices and T(a:, y, z, e) = (Tj;,Tj,,T^). Thus, the probabil¬ 
ity of finding odd-frequency triplet superconducting correla¬ 
tions with zero spin projection along the z-axis is 
Likewise, if T^; or Ty are finite, there exists triplet corre¬ 
lations with ±1 spin projections along the spin quantization 

35,36,82-84 

We have normalized all lengths by the superconducting co¬ 
herence length ^s- The quasiparticles’ energy e and the ex¬ 
change energy intensity are normalized by the zero tempera¬ 
ture superconducting order parameter, Aq. We also assume a 
low temperature of T = 0.05Tc, where Tc is the critical tem¬ 
perature of the bulk superconducting banks. We use natural 
units, with h = ks — 1, where fcs is the Boltzmann con¬ 
stant and define /3 = kBT/2. We consider weak exchange 
field strengths of \h\ = 5Ao corresponding to that found in 
ferromagnet alloys^^ such as, e.g., Cua;Nij^. A barrier resis¬ 
tance of C = 4.0 ensures sufficiently opaque S/F interfaces 
leading to appropriate solutions to the Usadel equations within 
the low proximity limit, F{x, y, z, e) <C 1. Having now out¬ 
lined the Keldysh-Usadel quasiclassical formalism and spin- 
parametrization framework employed in this work, we now 
proceed to present our analytical and numerical findings in 
the next sections. 


III. ONE-DIMENSIONAL HYBRID STRUCTURES 

In this section, we consider the effectively one-dimensional 
hybrid structures sketched in Fig. 1. We first derive the Us¬ 
adel equations and boundary conditions. We then decompose 
the supercurrent into its singlet and triplet components using 
the three-dimensional spin-parametrization scheme discussed 
in Sec. II.2. Utilizing this singlet-triplet decomposition, we 
analyze and characterize the supercurrent behavior based on 
the individual components involved. 


III.l. Spin-parameterized supercurrent 

Upon decomposing the Green’s function via Eq. (9), the 
Usadel equation, Eq. (1), transforms into the following 
eight coupled complex partial differential equations for one¬ 
dimensional systems, 


D e) -F e)} -F i {— 2£(q=Ta;(—e) -F iTy(—e)) 2S(— £)(/ix T — 0 (10^) 

D {Ta^S(-£) + dlTz{-e)) + i {t2T,(-£)L, t ‘FTy{.-£)hy - 2(TS(-e) + T,(-£))(£ ± /i.)} = 0 (10b) 

D {Ta^n(£) - idlTlie)] + i {±2{hx ± ihy)^*{e) - 2£(tT:(£) - iT;(£))} = 0 (10c) 

D {Ta^S*(£) + dlTlie)} + i {2{-e ± + T:(£)) ± 2hxT:{e) ± 2hyT;{e)} = 0 (lOd) 


The Kupriyanov-Lukichev boundary conditions at the left ing to the following differential equations: 

S/F interface, Eq. (4), are transformed in the same way, lead- ^ ^ 

iCdx - c*(£))(TS(-e) +T,(-£)) = Ts*(e), (Hb) 
- c*(£))(tT:(£) - 2 T;(£)) = 0, (11c) 

iCdx - c*(£))(tS*(£) + T:(£)) = ±S*(£). (lid) 
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Similarly, the Kupriyanov-Lukichev boundary conditions at 
the right F/S interface are also transformed as, 

iC^x + C*(£))(Tl'a;(~£) + —£)) = 0, (12a) 

(Ca. + c*(e))(TS(-£) + T,(-e)) = ±s*(e), (12b) 
(Ca. + c*(£))(TT:(e) - ^T;(e)) = 0, (12c) 

(Ca. + c*(e))(TS*(£) + T:(e)) = Ts*(e). (12d) 

By solving this coupled set of complex differential equations 
with the boundary conditions [Eqs. (11) and (12)], the rel¬ 
evant physical quantities can be obtained. We consider the 
cc-axis be normal to the interfaces, as shown in Fig. 5. The de¬ 
composition introduced above leads to the following expres¬ 
sion for the supercurrent density within the junction [Eq. (6)]: 

/ OO 

de { 

-OO 

S{e)d^S*{-e) - S(-e)5^S*(e) -f S{e)*d^S{-e)- 
Si-eyd^Sis) - dyTy-e)Tl{e) + dyT^{e)Tl{-e)- 
dy^l{-e)Tye) + d,Tl{e)Ty-e) - dy^y{-e)Tl{e)+ 
dyTy{e)Tl{-e) - d,Ty(-e)Ty{e) + 5,T;(e)T,{-e)- 
dyfy-e)Tl{e) + dy^ye)T*y-s) - d,T*y-e)Tye) + 
5a;T*(£)T2(-£)}tanh(e/3). (13) 


The current through the junction can be easily obtained by 
integration of the current density along the y and z directions 
over the F junction widths W and L, respectively (corre¬ 
sponding to the cross section of the wire). We assume that 
our system is very wide in the y direction, so that the one¬ 
dimensional approximation is valid, and therefore the current 
density remains constant in the y and z directions. It is con¬ 
venient to define the normalization constant, /g = LWJq, 
for the supercurrent, I{x). To extract the contributions to the 
total supercurrent from the even frequency singlet, and odd 
frequency triplet correlations, we have also decomposed the 
supercurrent accordingly into four components; 


Isoix) = lo / ci£{S(£)9,S*(-£) - S(-£)9,S*(£)+ 

S{ey dx'S{—e) — S(—£)*9a;S(£)} tanh(£/3), (14a) 

/ OO 

de {-dxTx(-e)T*{e) + dxTx{e)T*y-e) 

-OO 

-93;T*(-£)T;j(£) -f aa;T*(£)T;j(-£)}tanh(£/3), (14b) 

/ OO 

de {-5a;Ty(-£)T*(£) -f dxTy{e)T*{-e) 

-OO 

-a,T;(-£)Ty(£) +a,T;(£)Tj,(-£)}tanh(£/3), (14c) 

/ OO 

d£{9,T,(-£)T:(£) +a,T,(£)T:(-£) 

-OO 

-a^T*(-£)T^(£) -f i9a;T*(£)T2(-£)}tanh(£/3), (14d) 


where the total supercurrent is thus the sum of decomposed 
terms, namely. 


Itot{x) = Iso{x) + Isx{x) + Isy{x) + Isz{x). (15) 

This decomposition allows for pinpointing the exact behavior 
of the even- and odd-frequency supercurrent components. 


Results and discussions 

Various analytical or numerical schemes with varying ap¬ 
proximations have been employed to investigate the structures 
with magnetization patterns shown in Fig. 115,24,25,42^5 jjj 
analytical treatments"^^’*^, limiting approximations were em¬ 
ployed. For example, to study the noncollinear S/F/F/FIS 
structures^^45^ transparent boundaries are employed at the 
S/F interfaces together with the assumption that the anoma¬ 
lous Green’s function varies enough slowly through the mag¬ 
netic trilayer to warrant its Taylor expansion. Our full numeri¬ 
cal results involve no such approximations, and hence reveals 
cases where the inclusion of such effects may be important. 
One of the main aspects that our numerical approach reveals 
is the crucial role that each of the different types of super¬ 
conducting correlations play in the total supercurrent, given 
by Eq. (6). The supercurrent is composed of different com¬ 
ponents of even-frequency singlet and odd-frequency triplet 
correlations S, T, respectively. Such a decomposition is often 
neglected in Josephson structures that involve intricate mag¬ 
netic textures. 

For comparison purposes, we first consider the simpler 
S/F/F/S heterostructureshown in Fig. 1(a). Most of the 
results have dpi ^ dF 2 , and in all cases shown, the magne¬ 
tization of the right F layer is fixed along the z-axis while 
the left F layer magnetization rotates in the yz plane, that is, 
hi = ho(0, sin/3i, cos/?i). The magnetization orientation of 
the left F is thus characterized by the angle jdi, since the mag¬ 
netization is entirely in-plane. The top row of Fig. 2 illustrates 
the total critical supercurrent Itot and its decomposed com¬ 
ponents (/so, Isx, Isy, Isz) versus the thickness of right F 
layer, dF 2 - The current components generally vary with posi¬ 
tion X, so in order to display an overall view of their behavior 
as a function of magnetization orientation, we spatially aver¬ 
age [denoted by (...)] each component over dF 2 - The thick¬ 
ness of the left F layer is set typically at dpi = O.lS^s, and 
its magnetization has two components hy, and hz, using a rep¬ 
resentative angle of /3i = 0.27r (these values are chosen to in 
part support our comparison purposes in Subsec. 1V.6). To 
show the fine features of the O-tt transition profiles, we have 
used a logarithmic scale for the magnitude of the critical su¬ 
percurrent and its decomposed components. The critical cur¬ 
rent (far left panel) undergoes multiple O-tt transitions when 
varying dF 2 - The decomposed current components {Iso) and 
(Isz) also show the same behavior as seen in the remaining 
panels. Next, in the bottom row of Fig. 2, we plot the max¬ 
imum current and its components as a function of position 
for two representative values of the right F layer’s thickness 
idF 2 = 0-25,^s, and dF 2 = l-05^s). The current compo¬ 
nents, in contrast to the total current, often vary inside the 
magnetic layers: Iso, and Isz are shown to propagate within 
the two F layers. The spin-1 triplet component, Isy, is how¬ 
ever localized within the left F layer where /32 — 0.27r and 
thus hy y 0. We have investigated a wide range of parameter 
sets, involving /3i, and the superconducting phase differ¬ 
ences. The Isy component, does not propagate into the right 
F region where the exchange field is directed along z, demon¬ 
strating consistency with previous studies'^^. As discussed in 
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FIG. 2. (Color online) Top row. Critical supercurrent and its components as a function of dF 2 , the width of the right magnetic layer in a 
S/F/F/S junction (see Fig. 1(a)). The width of the first magnetic layer is fixed at dpi = O.lS^g. The current components are calculated 
by performing a spatial average over the right F layer / lo), {I Syix) / Iq) , {Isz{x) / Iq)). Bottom row: Spatial behavior of the critical 

Josephson current and its components as a function of position x inside the junction for two typical values of dF2 = 0.25^s, and d = l.OS^g. 
The magnetization orientation in the left F layer is fixed at a representative value of /?i = 0.27r, while the magnetization in the right F layer 
always points along z, corresponding to /32 = 0. 


the introduction, recent theoretical works showed that signa¬ 
tures of the triplet supercmi'ent may be detected by the ap¬ 
pearance of a second harmonic in the supercmi'ent in ballis¬ 
tic S/F/F/S Josephson junction, provided that dpi ^ dp 2 . 
The higher harmonics were shown to decay exponentially 
(faster than the first harmonic) when varying the system pa¬ 
rameters such as the thickness of the magnetic layers, and 
exchange field intensities, in the full proximity limit of the 
diffusive regime.^ Therefore, in the low proximity limit we 
consider in our manuscript, the higher harmonics are absent.^ 
It has been suggested that a trilayer"^^’"^^ of uniform magnetic 
materials with noncollinear magnetizations can reveal the sig¬ 
natures of long-ranged spin-triplet correlations where the two 
outer layers produce nonzero spin projections which can be 
detected in the middle F layer with orthogonal magnetiza¬ 
tion. To elucidate the source of the long-range triplet behavior 
in these types of trilayer configurations"^^, we investigate next 
the details of the individual components comprising the total 
supercurrent. 

We therefore consider a SIF/F/F/S trilayer structure, as 
depicted in Fig. 1(b). The magnetization of the central F 
layer is pinned in the z orientation, coinciding with the spin- 
quantization axis. The relative in-plane magnetization direc¬ 
tions in the sutTounding left and right F layers are described 
simply by the angles /3i, and j32, respectively. We denote the 
thicknesses of the left, middle, and right F layers by dpi, 
dp 2 , and dp^, respectively. In the top set of panels in Fig. 3, 
the spatially averaged supercurrent and its singlet and triplet 
constituents are shown as a function of dp 2 . We here show the 
results for dpi = dp^ = which is representative of the nu¬ 
merous equal-width cases investigated numerically. The bot¬ 
tom set of panels illustrate the spatial behavior of the super- 
curi'ent and its components. To isolate the spin-1 triplet con¬ 
tribution, Ipy, to the supercurrent in the middle F layer, we 
set j3i = Ps = it/ 2. In this case, the magnetization of the two 


outer F layers are strictly along y and orthogonal to the ex¬ 
change field direction of middle F. Such a magnetization con¬ 
figuration has been suggested as optimal for detecting the sig¬ 
natures of the spin triplet supercurrents"^^. As seen in the fig¬ 
ure, the critical supercurrent versus the middle F layer thick¬ 
ness shows multiple O-tt transitions, corresponding to points 
where the current nearly vanishes and then eventually changes 
sign. The averaged components of the total supercurrent. Iso 
and Isz demonstrate short-range signatures as exhibited by 
the multiple cusps compared to the equal-spin triplet compo¬ 
nent Isy The total critical current behavior is dominated by 
the triplet term, Isy, which undergoes fewer sign changes than 
the other components, and consequently fewer O-tt transitions, 
when changing dp2. Thus for the regime considered here, the 
supercurrent does not exhibit a very slow monotonic decay as 
a function of the central magnetic junction thickness, as re¬ 
flected in the absence of long-ranged behavior in Isy vs dp2. 

To further explore the behavior of the current throughout 
the junction, we next examine (bottom row. Fig. 3) the spa¬ 
tial dependence to the total supercurrent and its singlet and 
triplet components for representative values of unequal mid¬ 
dle F layer thicknesses, dp 2 — 0.55^s, and dp 2 — l.lS^g. 
We immediately observe from the left panel that as expected, 
the maximum total supercurrent is a constant in all parts of the 
junction, reflecting conservation of current there. The triplet 
component with zero-spin projection, Ipz, is localized in the 
middle F layer where the magnetization is directed along z. In 
contrast, the singlet component Iso oscillates throughout the 
junction while the triplet component, Isy, propagates without 
decay in the middle ferromagnet, which has its magnetiza¬ 
tion direction orthogonal to the spin orientation of Isy Thus, 
we observe that Isy is long-ranged in the middle F layer. 
Although Isy is spatially constant throughout the middle F 
layer, it changes sign, depending on dp2. This characteristic 
is seen in the top row of Fig. 3. We have investigated with our 
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FIG. 3. (Color online) Top row. Critical supercurrent and its components against the thickness of the middle F layer dp 2 in the S/F/F/F/S 
structure shown in Fig. 1(b). The thickness of the outer F layers, dpi and dps., are equal and fixed at ^s. The components are spatially 
averaged over the middle F layer thickness ({1sq{x)/1o), {lsy{x)/I o), {Isz{x)/I q}). Bottom row. Critical current as a function of position 
inside the junction x for two different thicknesses of the middle F layer, dp 2 = O.SS^s, and dp 2 — l.lS^s. The magnetization of the 
middle F layer is fixed along the z direction, P 2 = 0, while the magnetization is oriented towards the y direction in the outer F layers 
(/3i = /?3 = 7r/2). 


full numerical method several different geometrical parame¬ 
ter sets, including, e.g., much smaller dpi, and dps, as well 
as other jdi, and f}^. We typically found that the results pre¬ 
sented in the top row of Fig. 3 are quite representative of the 
singlet and triplet supercutTent behavior for the low proximity 
regime. 

In a recent experiment^® involving Ho/F/Ho trilayers 
with inhomogeneous magnetization patterns, a long-ranged 
Josephson supercurrent through the ferromagnet was detected. 
We here fully characterize the long-ranged triplet correlations 
in such S/Ho/F/Ho/S Josephson junctions (see Fig. 1(c)). 
The sandwiched central F layer represents a material with 
uniform magnetization, e.g., Cobalt, while the outer two Ho 
layers represent ferromagnets with a conical magnetization 
texture, such as that found in Holmium {Ho). As with the 
previous structures, we focus the study on the supercurrent 
behavior as a function of the central F layer thickness. We 
find that the supercurrent decays uniformly without any sign 
change when varying the middle F layer thickness, in agree¬ 
ment with other works"^^. However, the precise underlying 
role of the singlet and triplet components in the overall super- 
cutTent behavior has been lacking. In Fig. 4, we exhibit the 
total current and its decomposition for a variety of parameters 
appropriate for the inhomogeneous conical magnetic junction. 
The middle F layer is magnetized along z, while the con¬ 
ical magnetization patterns of the Holmium layers have the 
adopted form: 

hpo = ho (cosa,sinasin( 7 x/a),sinacos( 7 x/a)), (16) 

with the following material parameters: a = 0 . 02 (i is the dis¬ 
tance of interatomic layers, d denotes the junction thickness, 
a = 4 tt/9 is the apex angle, and 7 = 7r/6 is the rotation angle 


of the cone structure, consistent with experimental values^^. 
The magnitude of the exchange field is unchanged through¬ 
out the F regions. The top set of panels in Fig. 4 shows the 
critical total current and its components as a function of the 
middle uniform F layer thickness dp. The components Iso, 
Isy, and Isz are averaged over the middle F layer [denoted by 
(...)]. We assume the two Ho layers have identical magnetiza¬ 
tion patterns"^^, and both thicknesses equal dpo = The 

first panel on the left shows the total critical supercurrent ver¬ 
sus dp and exhibits the expected decay over a few coherence 
lengths. Examining the signatures of the other components in 
the top panels, we see that the opposite-spin singlet and triplet 
components. Iso and Isz respectively, demonstrate well de¬ 
fined oscillatory behavior and corresponding sign changes as 
a function of dp. We can also conclude that the net super- 
current arises mainly from the spin -1 projection of the triplet 
components, Isy, which is long-ranged in the middle F with 
an exchange field direction orthogonal to the spin-orientation 
of Isy. For the Holmium magnetization profile, hx does not 
vary in space, and therefore Isx behaves similarly to the Isy 
component in the S/F/F/F/S structure (see Fig. 3). Al¬ 
though Isx undergoes fewer sign changes when varying dp, 
the O-TT transitions are clearly present for this component. The 
triplet component, Isy, the main contributor to the total cur¬ 
rent, does not switch directions when increasing the middle 
F layer thickness dp, and its magnitude often dominates the 
other components. Its nearly monotonic decay can be traced 
back to the corresponding component of the magnetization 
profile: in the Ho layer, hy rotates sinusoidally as a func¬ 
tion of position, generating long-ranged odd-frequency cor¬ 
relations that are not subject to the spin-splitting effects of 
the magnet responsible also for the oscillatory behavior of su- 
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FIG. 4. (Color online) Top row. Critical current and its components versus middle F layer thickness, dp in S/HojFjHo/S structure 
where the outer F layers’ magnetization pattern is Holmium (ii’o)-like . The configuration is illustrated schematically in Fig. 1(c). The 
thicknesses of the outer F layers are fixed at dpo = 1.4^s, whereas the magnetization is fixed along the 2 ; direction in the middle F 
layer, j52 = 0. The components of the critical supercurrents are calculated hy spatially averaging over the thickness of the middle F layer 
({Iso{x)/Io), {Isy{x)/Io}, {Isz{x)/Io)). Bottom row: Critical supercurrent against position inside the magnetic junction x where the thick¬ 
ness of middle F layer has two values dp = 0.45^s, 1.2^s. 


perconducting correlations. We have also found consistency 
with previous studies"^^, where the monotonic decay of the su¬ 
percurrent appears when dno is large enough to contain at 
least one spiral period. On the other hand, the sign chang¬ 
ing behavior emerges for small dno, so that the Ho layers 
effectively mimics a uniform ferromagnet. This aspect was 
revisited in a recent work employing a lattice model.*® The 
bottom row of Fig. 4 illustrates the decomposed components 
of the maximum total supercurrent as a function of position 
throughout the ferromagnet regions. Two representative thick¬ 
nesses of the middle F layer are considered; dp = 0.4^5, and 
dp = 1 . 2 ^ 5 . Similar to the S'/T’/F/F’/S'junction above, the 
triplet components with spin projection m = ±1 on the z-axis 
{Isx and Isy), are constant over the entire middle F region 
and therefore can be classified as long-ranged. To summarize 
this section, we studied the behavior of the critical supercur¬ 
rent through low-proximity one-dimensional structures shown 
in Fig. 1. By directly decomposing the supercurrent using the 
spin-parametrization technique given in the theoretical meth¬ 
ods section, we numerically studied the origins of the super¬ 
current behavior in terms of its short-ranged and long-ranged 
components. Our results showed that S/F/FIS structures do 
not support any long-ranged supercurrent components, reaf¬ 
firming the findings of Ref. 45, while S/F/F/F/S junctions 
host long-ranged supercurrent components'^^ that are more 
prominent in the inhomogeneous SIHolFIHolS structures 
as experimentally observed in Ref. 50 and verified theoret¬ 
ically by numerical studied of full proximity regime in Ref. 
43. We showed that the long-ranged supercurrent component 
corresponds to the rotating component of the magnetization 
texture, and is the main contributor to the total supercurrent. 
The numerical results presented in this section shall be used 


for later comparisons when we categorize structures into two 
classes based on the supercurrent direction with respect to the 
FIF interface orientation; parallel or perpendicular. The one¬ 
dimensional structures in Fig. 1 belong to the latter class. We 
now direct our attention to two-dimensional hybrids, includ¬ 
ing the possibility of an applied magnetic field. The singlet- 
triplet decompositions discussed above shall be employed to 
pinpoint exactly the spatial behavior of the associated compo¬ 
nents of the total charge supercurrent. 


IV. TWO-DIMENSIONAL HYBRID STRUCTURES 

In this section we present the main results of the paper. 
We first consider a two-dimensional magnetic SIF/S system 
subject to an external magnetic field for two regimes; the 
wide junction Wp ^ dp, and the narrow junction Wp <C dp 
regimes. The magnetic strips are sandwiched between two s- 
wave superconducting reservoirs, where the exchange field in 
the F strips is spatially uniform. This permits analytical solu¬ 
tions that are capable of accurately predicting the behavior of 
relevant physical quantities such as charge and spin supercur¬ 
rents, as well as the pair potential. The analytical results are 
complemented with full numerical investigations, including 
studies of the dependence of the critical charge supercurrent 
on the external magnetic field, and the corresponding appear¬ 
ance of proximity vortices. The current density spatial map 
is also investigated, giving a global view of the distribution 
of supercurrents throughout the junction. We next consider 
two kinds of finite-sized magnetic SIF/F/S Josephson junc¬ 
tions subject to an external magnetic field. For these systems, 
analytical routes are scarce, and we must in general resort to 
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numerical approaches. In one case, we assume the double 
layer magnetic F strips comprising the F/F junction are par¬ 
allel with the S interfaces. In the other case, however, we 
assume that the F strips are perpendicular to the S interfaces. 
Using the spin-parametrization introduced in Sec. II.2, we can 
then study the even and odd frequency components of the total 
charge Josephson current inside the proposed structures. 

IV. 1. Technical assumptions and parameters 

In this subsection, we discuss the assumptions used in 
our calculations, along with the parameters and notations 
used throughout. As was previously mentioned, the exter¬ 
nal magnetic field is confined within the non-superconducting 
regions®^’’(see Fig. 5). We restrict the magnetic field to 
be oriented perpendicular to the junction plane, which for 
our coordinate system corresponds to the z-axis. Thus, su¬ 
percurrent flow resides in the xy plane. The situation where 
the external magnetic held is parallel to the current direc¬ 
tion has been studied both experimentally and theoretically. 
We can therefore represent the magnetic held as, H{x) = 
HQ{x)Q{dp — x), where dp is the junction length and H = 
(0, 0, Hz). This assumption also ensures that the macroscopic 
phases of the two superconducting electrodes are unaffected 
by the external magnetic held. This widely used assumption 
has demonstrated good qualitative agreement with experimen¬ 
tal measure ments.®*’^*’^^’^® The Josephson junctions investi¬ 
gated in this work therefore are assumed to have negligible 
magnetic held screening.If on the contrary, the magnetic 
held is not restricted to the F regions, it becomes necessary 
to solve a set of partial differential equations, Eq. (1), self- 
consistently in tandem with Maxwell’s equations and the su¬ 
perconducting order parameter A(a:, y, T). Hence, a suitable 
choice for the vector potential A that we use satisfying the 
Coulomb gauge, V • A = 0, is, A{x) = {—yHz{x),Q,A). 



FIG. 5. Proposed setup of the Josephson S/F/S junction subject to 
an external magnetic field H. The magnetic wire (F) with rectan¬ 
gular dimensions, dp and Wp, is sandwiched between two s-wave 
superconducting (S) leads. The exchange field of the magnetic layer, 
h = (h^, , h^), can take arbitrary orientations. The external mag¬ 

netic field is oriented along the z direction. The two-dimensional 
system resides in the xy plane so that the interfaces are along the y 
direction and x axis is normal to the junction. 


In normalizing our equations, we write the external magnetic 
flux as $ = where $ = WpdpHz, Wp is the junction 

width, and $o = ^/2e is the magnetic flux quantum. 


IV.2. Uniform S/F/S heterostructures 

We here consider a two-dimensional fetTomagnetic S/F/S 
Josephson junction where the magnetization of the magnetic 
strip is homogeneous. Although our theoretical approach al¬ 
lows for completely general patterns in the magnetization tex¬ 
ture, we restrict our focus here to a specific case where the 
exchange held has only one component along the z direction, 
h{x, y, z) = (0, 0, h^), thus permitting analytical solutions to 
the Usadel equation. The two-dimensional junction resides in 
the xy plane so that the S/F interfaces are parallel with the y 
axis (see Fig. 5). The cotTesponding system of coupled partial 
differential equations [Eq. (1)], now reduces to a smaller set 
of decoupled partial differential equations. We are then able 
to derive analytical expressions for the anomalous component 
of the Green’s function and therefore the charge supercurrent 
and pair potential. 

Using this simplified system of decoupled partial differen¬ 
tial equations, we consider two regimes: In the first case, we 
assume the junction width Wp ^ dp, thus, terms involving 
the ratio dp/Wp can be dropped, leading to further simplifi¬ 
cations. In the second regime, the junction width Wp dp, 
corresponding to a narrow magnetic nanowire. To be com¬ 
plete, we also implement a full numerical investigation, with¬ 
out the simplifying assumptions above and with arbitrary val¬ 
ues of ratio Wp/dp. This requires numerical solutions to a 
complex system of partial differential equations [see Eqs. (1)]. 
Several checks on the numerics were performed, including 
reproducing previous results involving nonmagnetic S/N/S 
Josephson junctions, where the exchange held of the F layer 
is equal to zero.®^’^* 

The full Usadel equations in the presence of an external 
magnetic held H, and corresponding vector potential. A, are 
written: 

VV±(-e) - 2ieV • A/±(-e) - AieA ■ V/±(-e) - 
Ae^A^Ui-e) - /«(-£) = 0, (17a) 

+ 2*eV • If lie) + 4*eA • VfHe) - 
4e^A^flie)-^^^^^^flie) = 0, (17b) 

where V = {dx, dy,dz). The above decoupled partial differ¬ 
ential equations appear only for a magnetic junction where the 
magnetization has one component h = (0,0, h^). If we now 
expand the boundary conditions given by Eq. (4), at the left 
S/F interface, we And, 

{C(V - 2ieA) - c*ie) ± = Ts*(e)e-(!'§a) 

{C(V + 2ieA) - c*ie) T = ±s*(e)e*'',(18b) 
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while the boundary conditions at the right interface take the 
following form: 

{C(V - 2ieA) + c*{e) T = ±s*(e)e-t^’‘^a) 

{C(V + 2ieA) + c*{e) ± i^}/;(e) = Ts*(e)e*®^.(19b) 

The simplifying geometric approximations mentioned above 
can now be applied to the above equations, while adhering 
to the requirement that the corresponding regimes are experi¬ 
mentally accessible. 


Wide junction limit, Wf dpt Analytical results 

If we assume that the width of junction is much larger than 
its length, terms involving dp/Wp in Eq. (17) can be ne¬ 
glected, yielding the following decoupled Usadel equations: 


-2i"^/±(-e) =0, (20a) 

Ct 

'9^/±(e) + 4i$y5^/±(e) - 

-2i^—fli£)=0. (20b) 

ep 

Here, we define e± = e ± h^{x,y,z), and the Thouless 
energy ep = D/d^. Note that all partial derivatives in 
Eqs. (20) are solely with respect to the x coordinate. In other 
words, the original two-dimensional problem is now reduced 
to a quasi one-dimensional one. These uncoupled differential 
equations can be solved analytically thus permitting additional 
insight into the transport properties of ferromagnetic Joseph- 
son junctions. The Kupriyanov-Lukichev boundary condi¬ 
tions, Eq. (4), at the left StF interface located at a; = 0 re¬ 
duces to: 

|c(5x + 2iy$) -f c(-e) T *^|/±(-£) = Ts(-£)d^'l4) 

|C(9x - -f c(-e) ± i^}/i(e) = ±s(-£)e"(flb) 

Similarly, the boundary conditions at the right SIF interface 
located at x = dp can be written as follows: 

|C('9a; + 2%$) - c{-£) ± i^|/±(-£) = ±s(-e)et2^q) 

{C(9, - 2%$) _ c(-£) T *^}/i(e) = Ts(-e)e*^22b) 

The macroscopic phases of the left and right superconducting 
terminals are labeled di and dr, respectively. The magnetic 
strips are assumed isolated in the y direction so that physi¬ 
cally no cutTent passes through the boundaries at y = 0, and 
y = Wp. Thus, to ensure that the supercurrent does not pass 
through the vacuum boundaries in the y direction, we have the 
following conditions: 

(23a) 
(23b) 


With the solutions to Eq. (20) at hand, we are now in a po¬ 
sition to calculate the current density and the pair potential for 
a given magnetic flux. The current density, given by Eq. (6), 
can thus be expressed as: 

/ OO 

d£ tanh(e 

-OO 

+f+{-£)Vr_{£) - f+{£)Vr_{-£) - /_(£) 

v/;(-£) + f*_{-£)vf+{£) + /;(-£)v/_(e) 

-/;(e)V/_(-e) - f*_{e)Vf+{-£) + UeA[ 


/3){/_(-e)V/;(e) 


r (£)/+(-£) - /;(-£)/-(£) - fl{-s)U{£) 

+ /;(£)/_(-£)]}. (24) 

The normalized pair potential, Eq. (8) now reads: 

/ OO 

[f+i^) - /-!-(-£) + 

-OO 

/- (-e) - /- (e)] tanh(£/3)d£. (25) 

If we solve the Usadel equations Eqs. (20) using the boundary 
conditions (21) and (23), we arrive at the following solutions 
to the anomalous component of the Green’s function: 

f±i-£) = (26) 

where the numerators and denominators are given by 

= s*(£)e-5*(^"^'^+^){e"*^^[CA^cosh(xA^) -f 

(c*(£) ± sinh(a;A;F)] -|-e*‘^[CA=p x 

LtJ' 

cosh(Azp — ccAzp) -I- (c*(£) ± sinh(A;p — xAzp)] |, 


= sinh(Azp)|c^A^ -f (c*(£) ± + 

2CAzp(c*(£) ±i^) cosh(Azp). 

Similar solutions can be found for /^(—£). In order to sim¬ 
plify notation, we have defined A± = 2z(e ± h^)/ep, and 
the macroscopic phase difference of the superconducting ter¬ 
minals is denoted hy (j> = 6i — dr- The charge cutTent den¬ 
sity expressed by Eq. (24) involves eight terms, f± (±£), and 
/^(±£), which should be derived to obtain an analytical ex¬ 
pression for the supercurrent flow. In our calculations thus far, 
the interfaces are assumed spin-active, namely, Gp/Gp ^ 0. 
To maintain tractable analytic solutions, we drop the c(£)- 
terms in addition to the spin-active contributions, which is 
appropriate for experimental conditions involving highly im¬ 
pure superconducting terminals. These widely used approx¬ 
imations lead to intuitive and physically relevant analytical 
solutions. ^4.24 Substituting the solutions to the Usadel equa¬ 
tions into Eq. (24), we find the charge supercurrent density in 
the X direction: 


m = Jo / d£ 


2itanh(£^)N 


- 1/2 


dysm{(j> — 2y$), 
(27) 


dyf±{±£) = 0 , 

dym±£) = o. 


— OO 


C'a_a+ 



12 


where, 

N = s*^(—e)^A+ csc(A_) + A_ csc(A+)^ + 

s*^(e) ^A+csch(A_) + A_csch(A+)^. 

Integrating the junction width over the j/-direction, we end up 
with the total charge supercurrent across the junction; 

sin $ 

/((/),$)=/o sin (28) 

where we have extracted the phase and flux dependent terms 
and absorbed the remaining coefficients into Ig. The max¬ 
imum charge supercurrent occurs when the superconducting 
phase difference equals \4)\ = 7r/2. From Eq. (28), we see 
immediately that the critical charge current exhibits the well- 
known Fraunhofer interference diffraction pattern as a func¬ 
tion of the externally applied flux $. We also recover the re¬ 
sults of a normal 5'/-/V/S'junction®*, where = 0. Thus, our 
analytical expressions for wide S/F/S Josephson junctions 
experiencing perpendicularly directed external magnetic flux 
yields the same critical current response as a normal S/N/S 
junction®*). Here, however, there are additional, experimen¬ 
tally tunable physical quantities which can cause sign changes 
in Ig{h^ ,G s/G t), and consequently J(^, fl)). When present¬ 
ing a global view of the current density, it is illustrative to 
examine a spatial map of its behavior. By utilizing Eq. (6), it 
is possible to calculate the charge current density throughout 
the xy plane for a wide magnetic S/F/S Josephson junction. 

If we now insert the recently obtained solutions to the Us- 
adel equations into the pair potential equation, Eq. (25), we 
arrive at the following analytical formula which provides a 
spatial map of [/pah for a wide junction; 

/ -too 

del e“ 5([s*(—e)A+ cos(a;A_) x 

-OO ^ 

csc(A_) -I- s*(—e)A_ cos(a;A+) csc(A_|_) -I- s*(e)A+ x 
cosh(xA_)csch(A_) -I- s*(£)A+ cosh(a;A_|_)csch(A+)]e^**'* 


}i^ = 0, recovers the nonmagnetic S/N/S junction result for 
the proximity pair potential.®* Comparing with the current 
density, Eq. (27), we see that the current density and prox¬ 
imity pair potential both vanish at the same locations, how¬ 
ever the current density vanishes at additional positions cor¬ 
responding to sin — 2?/$) = 0, or, when (jj — 2j/$ = rmr 
for m = 0, ±1, ±2, ±3... . The origin of these extra zeroes 
in the current density arises from the cancellation of counter- 
propagating currents from the orbital motion of the quasipar¬ 
ticles. These paths are visualized using spatial mappings, pre¬ 
sented below. 


Narrow junction limit, Wf ^ dp-' Analytical results 

The next useful regime that leads to analytical results is 
that corresponding to a narrow junction, that is, Wp dp. 
In this case, we assume that the width of the ferromagnetic 
layer Wp < ^h, where is the characteristic length de¬ 
scribing the Green’s function oscillations in the ferromagnetic 
layer. This assumption permits averaging the relevant equa¬ 
tions over the junction width. Therefore, making the substi¬ 
tutions {—Hyx)y = 0 and {F[^y^)y = H^/12, where (.. .)y 
denotes spatial averaging over the y direction, results in the 
modified Usadel equations; 

dlf±{-e) - = 0, (31a) 

Ct 

dlfUe)-2fi^^^)fUs) = 0, (31b) 

ep 

where we define Tp = Dn^Fl^Wp/G, and thus, Vpjtp = 
$^/6. The quantity, E/f, is known as the magnetic depairing 
energy®*. As can be seen, the above assumptions have consid¬ 
erably simplified the Usadel equations, and we are now able 
to efficiently solve the differential equations and derive ana¬ 
lytical expressions for useful physical quantities. After some 
straightforward calculations, we arrive at the following solu¬ 
tions to the Usadel equations for the anomalous Green’s func¬ 
tion; 


-l-[s*(— £)A+ cos(A_ — xA_) csc(A_) -I- s*{—e)X- x 
cos(A+ — a;A+) csc(A+) -I- s*(£)A+ cosh(A_ — xX-) x 

csch(A_) -I-s*(£)A_ cosh(A+ — a;A+)csch(A+)]e*‘^)| tanh(^^) /±(—e) = (32) 


If we restrict the proximity pair potential profile above by con¬ 
sidering a fixed x position, corresponding to the middle of the 
junction (a; = 1 /2), we arrive at an expression which is now 
X independent; 


[[pair — 


^^cos {[(f) - 2y$]/2) J- 


|s*(-£)A+csc(A_/2) 

-I- s*(—£)A_ csc(A+/2) -I- s*(£)A+csch(A_/2) 

-I- s*(£)A_csch(A+/2)| tanh(£/3). (30) 


where the numerators and denominators T)^ are now de¬ 
fined by the following expressions; 

= s*{e)e~^ |(c*(e) ± ^ 

(sinh(xA)j)) + sinh(A)[) — a:A)j))) -I- x 

cosh(A)j) — a;A)[)) -I- CA){) cosh(a;A)[))|, 

and. 


As seen, the cos ([(/> —2y$]/2) term is zero at </)—2?/$ = tott, 
for m an odd integer. Therefore, the zeros of the proxim¬ 
ity pair potential at the middle of magnetic strip are located 
at y = {(j) — m7r)/2$ so that —1/2 < y < 1/2. Setting 


= sinh(A|)|c^A|^ -f (c*(£) ± + 

2(X+{c*{e)±i^) cosh(A|). 
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Here, we have introduced a simplified notation, where we de¬ 
fine; At^ = and _ jf substitute 

± ex’ ^ €x 

these solutions to the anomalous Green’s function, Eq. (32), 
into the Josephson current relation, Eq. (24), we arrive at the 
flow of charge supercurrent through a narrow S/F/S junction 
in the presence of an external magnetic flux $: 


m=j_ 


~ , 2itanh(£/3)N . , 
de— - ^—r sin(0), 

-oo C2AIA;A^A+ 


(33) 


N = —Al|s*(—£)^AlA+csch(A;j]) — 
s*(e)^Aj] (A;j]csch(At) -f Alcsch(A^)) | 
s*(—e)^Aj^AtA^csch(Al), 


where the junction characteristics are obtained via averaged 
values. The critical supercurrent thus shows a monotonic de¬ 
caying behavior against the external magnetic flux T* in the 
narrow junction regime. This fact can be understood by not¬ 
ing the role of T// in the Usadel equations. In the narrow 
junction limit, the external magnetic field breaks the coher¬ 
ence of Cooper pairs and the Eraunhofer diffraction patterns 
for the critical current in the wide junction limit turns to a 
monotonic decay. The proximity vortices, which are closely 
linked to the Eraunhofer patterns, vanish in this regime due to 
the narrow size of junction width. This finding is also in agree¬ 
ment with the nonmagnetic SIN IS junction counterpart®^’^'. 
If we substitute the solutions from Eq. (32) into the pair po¬ 
tential, Eq. (25), we find that the zeros which appeared in the 
wide junction limit have now vanished in the narrow junction 
regime. We now proceed to compliment our simplified ana¬ 
lytical results with a more complete numerical investigation. 


Interference patterns, proximity vortices, and current densities: 
Numerical investigations 


Previously we utilized various approximations to simplify 
situations and permit explicit analytical solutions. To inves¬ 
tigate the validity of our analytic results, we solve numeri¬ 
cally the Usadel equations, Eq. (17) with appropriate bound¬ 
ary conditions found in Eqs. (18) and (19). We now retain the 
c{e) terms, and include a broader range of junction widths, not 
only those that are very narrow or very wide, but also interme¬ 
diate widths that are not amenable to an analytical treatment. 
We shall present results for the current density and for spa¬ 
tial maps of the pair potential in the junction subjected to an 
external magnetic field. 

Eigure 6 exhibits the results of our numerical studies for the 
critical supercurrent response to an external magnetic flux, $, 
in a SIFIS junction with a variety of ratios Wp/dp- The 
junction length and exchange field intensity are fixed at par¬ 
ticular values of dp = 2.0^s and \h\ = S.OAq respectively. 
The interface resistance parameter is set to C = 4.0. These 
choices ensure the validity of the Green’s functions in the low 
proximity limit; F{x,y,z,e) 0.1F_'^°{x,y, z,e). To more 
clearly see the effects in the scaled plots, we normalize the 
critical charge supercurrent by the maximum of this quan¬ 
tity max(/c) for each case separately. The figure clearly re¬ 
veals a full Fraunhofer pattern for the case of a wide junction 
Wp = 10^5. This pattern is consistent with the analytic ex¬ 
pression for the critical current given by Eq. (28). The critical 
supercurrent decays with increasing magnetic field and un¬ 
dergoes a series of cusps, indicating that Ic changes sign. The 
observed sign-change reflects aspects of the orbital motions 
of the quasiparticles, while the decaying behavior is indicative 
of the pair breaking nature of an external magnetic field as we 




FIG. 6. Magnetic interference pattern of an SIFIS Josephson junc¬ 
tion (see Fig. 5). The critical supercurrent /c('fi), is shown normal¬ 
ized by its maximum value, max(7c), and plotted against the exter¬ 
nal magnetic flux 4>. Three values of the magnetic layer width Wp 
are considered; lO^s, 2.0^s, and 0.2^s, whereas the length of F is 
assumed fixed at a representative value of dp = 2.0^s. 


FIG. 7. Spatial maps of the pair potential and current density in the 
presence of an external magnetic field applied normally to the plane 
of an SIFIS Josephson junction (see Fig. 5). We assume the junction 
is relatively wide, Wp = lO^s, while its length is sufficiently short 
to allow for current transport across the x axis, namely, dp = 2 . 05 s. 
The external magnetic flux is 4> = 34>o, and the macroscopic phase 
difference is corresponds to \4)\ = n/2. 
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explained previously. A relatively narrower junction width of 
Wp = 2.0^s is also investigated. As shown, the ideal Fraun¬ 
hofer pattern is now modified, yet still retains its trademark 
signature. This diffraction pattern transitions to a uniformly 
decaying behavior for the sufficiently narrow junction with 
width Wp = 0.2^s. This narrowest junction confines the or¬ 
bital motion of the quasiparticles and thus the pair-breaking 
by the external field plays a dominant role in the critical su¬ 
percurrent response. 

In Fig. 7 we plot the 2D spatial dependence of the pair po¬ 
tential and charge current density for a wide junction, where 
the larger width can more effectively demonstrate the orbital 
paths as they relate to the pair correlations and supercurrent 
response. To be consistent with Fig. 6, we set $ = S^hg, 
Wp = lO^s, and \(j)\ = 7r/2. As shown, three zeros in the 
pair potential appear at the middle of the junction (at a; = 0) 
as was found analytically in Sec. IV.2. The intensity of the 
external magnetic flux fl) determines the number of zeros and 
distance between neighboring zeros. The panel on the right 
corresponds to a spatial and vector map of the charge cur¬ 
rent density, revealing the circulating paths of the quasipar¬ 
ticles. If we set \4)\ = 7r/2 in Eq. (30), with the constraint 
— 1/2 < y < 1/2, we find that y = —1/12, 3/12, 5/12 cor¬ 
responds precisely to the zeroes of the pair potential found 
from the general numerical treatment, shown in Fig. 7. Like¬ 
wise, y = ±1/12, y = ±3/12, and y — ±5/12 gives the 
zeroes of the current density. As seen, the additional zeros in 
the current density corresponds to locations where the trajec¬ 
tories of two opposing circular paths overlap, thus canceling 
one another. The other y values giving zeroes, which corre¬ 
late with the same locations as the pair potential, are at the 
cores of the circulating paths. Such behavior is reminiscent of 
Abrokosov vortices, where the supercurrent circulates around 


normal state cores. The zeros in the pair potential may thus 
be viewed as proximity vortices.Abrikosov vortices, 
which carry a single magnetic quantum flux $o, are however 
intrinsic to type II superconductors.One of the criteria 
for categorizing the superconducting state of a material is the 
nature of these intrinsic vortices.®* Nonetheless, such proxim¬ 
ity vortices are generally geometry-dependent and rely on the 
mutual interaction between a magnetic field with the super¬ 
conductor, in contrast to intrinsic Abrikosov vortices.®^ 

IV.3. Spin-parametrization 

The Usadel equations in this section have dealt solely with 
the even-frequency superconducting correlations with spin- 
zero projection along the spin-quantization axis. This is 
due to the fact that we have only considered ferromagnetic 
strips with a uniform magnetization texture. For inhomo¬ 
geneous magnetization textures, the complex partial differ¬ 
ential equations become coupled and increase in number to 
eight in the low proximity limit. This number is doubled if 
the full proximity limit^^ is considered. Fortunately, it has 
been well established that the low proximity limit is suffi¬ 
cient to capture the essential physical properties of proximity 
systems such as the ones proposed in this paper. To inves¬ 
tigate the behavior of even- and odd-frequency correlations, 
we employ a spin-parametrization technique^®’^^’®^ that has 
been frequently used to study the characteristics of magnetic 

SyStemS.3-4.33.36, 44.45.49,53,83,92 

If we now substitute this decomposition of the anomalous 
Green’s function into the Usadel equation, Eq. (1), and con¬ 
sider a two-dimensional system, we end up with the following 
coupled set of differential equations in the presence of an ex¬ 
ternal flux $; 


D 


I T dlT,{-e) T ± idlTy{-e) + + 4z$j/(T9.T,(-e) ± id,Ty{-e)) 

± iTy(-e))| ± idj^^Ty - 2e(TT,,(-e) ± iTj,(-e)) t 2S(-e)(/i" ± = 0, 


(34a) 


d\ t dlSi-e) T 




\WpJ y 


dl^{-e) + dlT,{-e) + 




\WpJ y 


dlT,{-e) ± 4i$y(Ta,S(-e) ± a,T,(-e)) 


- 4 $ 2 y 2 (^g(_£) + T,(-e))| ± idl[ ± 2T,(-e)/i" t 2Ty{-e)hy - 2(±S(-e) ± T,(-e))(£ ±h^)}= 0, (34b) 

d| t dlTlis) T {^YdlTlie) - " 4*$y(Ta,T:(£) - zd,T;{s)) 


- - *T;(e))| ± id%[ ± 2(/i" ± z/i«)S*(e) - 2e{TTUe) - 2T;(e))} = 0, 

T dlS*{e) T [^ydp*{e) + ± (^^ydXi^) - 4z$y(Ta,S*(e) ± a,T:(e)) 

- 4 $ 2 j/ 2 (^g*(^) T*(e)) j ± i4|2(-e ± h^){T^*{e) + T*(e)) ± 2/z^T*(e) ± 2OT;{£)} = 0. 


(34c) 


(34d) 
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The spin parameterized boundary conditions at the left S/F interface, Eq. (4), now reads: 

{c(a, + 2i^y) - c*(e)} ( T T,(-e) + iT^(-£)) = 0, (35a) 

{C(a, + 2i<i>y) - c*(e) ± *^} ( T S(-e) + T,(-£)) ± s*{e)F^‘ = 0, (35b) 

{C(a, - 2i^y) - c*(£)} ( T T:(£) - zT;(£)) = 0, (35c) 

{C(a,-2i<i>2/)-c*(£)Ti^}(T §*(£)+T:(£)) Ts*(£)e-*®' =0, (35d) 

and for the right S/F interface, we have: 

{C(a, + 2^$2/) + c*(£)}(TT,(-£) + ^TJ,(-£)) =0, (36a) 

{C(a. + 2*4>2/) + c*(£)T*^}(TS(-£)+T,(-£)) Ts*(£)e*"’' =0, (36b) 

{C(a,-2z<i>2/) + c*(£)}(TT:(£)-zT;(£)) =0, (36c) 

{C(a.-2z<i>j/) + c*(£)±i^}(TS*(£)+T:(£)) ±s*(£)e-*®’- =0. (36d) 


Using Eq. (6), the spin-parameterized current density through the junction (along the x direction) in the presence of an external 
magnetic field is: 


Ux,y) =Jo J d£|s(£)a,s*(-£) - S(-£)9,S*(£) + S(£)*a,S(-£) - S(-£)*9,S(£) - 5 ,t,(-£)t:(£) + 9,t,(£)t:(-£) 

- 9,T:(-£)T,(£) + c),T:(£)T,(-£) - a,T,(-£)T;(£) + a,T,(£)T;(-£) - 5,T;(-£)T,(£) + 9,T;(£)T,(-£) 

- 9,T,(-£)T:(£) + 9,T,(£)T:(-£) - a,T:(-£)T,(£) + a,T:(£)T,(-£) - 4i$y [s(£)S*(-£) - S(-£)S*(£) 

+T.(-£)T:(£) - T,(£)T:(-£) + T,(-£)T;(£) - T,(£)T;(-£) + T,(-£)T:(£) - T,(£)T:(-£)] I tanh(£/3). (37) 


The y component of the supercurrent density, Jy{x, y), takes the same form as Jx{x, y), with the partial-a; derivatives replaced 
by derivatives with respect to y, and also the inclusion of the appropriate component of the vector potential. 

Incorporating the spin-parametrization, Eq. (9), into the expression for the pair potential, Eq. (8), we end up with the following 
compact equation: 


Unair 


J d£|s(£) — S(—£)| tanh(£/3). 


(38) 



FIG. 8. Proposed setup of an S/F/F/S Josephson junction where the double layer F region has its interfaces parallel with the S/F interfaces. 
The ferromagnetic layers and s-wave superconducting terminals are labeled by Fi, F 2 , and S, respectively. Our theoretical technique allows 
consideration of exchange fields with arbitrary orientations. We thus denote the magnetization of each F layer by hi — (/if, /if, /if), for 
i = 1,2. The width of the rectangular F strips are equal and denoted by ITV, while the length of each magnetic layer can differ and are 
denoted dpi, and dF 2 - The junction is located in the xy plane and all interfaces are along the y axis. An external magnetic field H points 
along the z axis, normal to the junction plane. 
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FIG. 9. Spatial maps of the total critical current density, pair potential, and supercurrent components in an S/F/FIS Josephson junction for 
the configuration depicted in Fig. 8. The supercurrent is decomposed into its even y)) and odd-frequency (J"^y{x, y), J^^y{x, y)) 

components. The two ferromagnetic layers have unequal lengths: dpi = 0 . 2 ^ 3 , and dF 2 = 1.8^s, while their widths are equal: Wfi = 
Wf 2 = 2.0^s. The vertical lines identify the two ferromagnetic regions labeled Fi and F 2 . In part (a), there is no applied magnetic field, 
while in part (b) the external magnetic field corresponds to a flux of = 3T>o. The critical current components in both cases, (a) and (b), are 
plotted as a function of 2 ;-position at four different locations along the junction width: y = 0.25^s, O.SO^s, 0.75^s, l.O^s. 


Thus, the pair potential involves only even-frequency compo¬ 
nents of the parameterized Green’s function, consistent with 
the presence of a singlet order parameter in the left and right 
superconducting banks. In the next subsections, we solve 
the above equations for two different magnetization config¬ 
urations, with and without an external magnetic field. The 
spin-parametrization outlined here will then delineate the var¬ 
ious contributions the singlet and triplet correlations make to 
the supercurrent. 


IV.4. S/F/FIS heterostructures: parallel F/F and S/F 
interfaces 

Employing the spin-parametrization with the Usadel equa¬ 
tions together with the boundary conditions given in the pre¬ 
vious section, we numerically study the proximity induced 
triplet correlations in a two-dimensional magnetic S/F/F/S 
Josephson junction. Here the F/F interface is parallel with 
the outer S interfaces. The configuration is depicted in Fig. 8. 
The system is subject to an external magnetic field H and 
we study the behavior of the corresponding spin-triplet cor¬ 
relations. We also compare these results with the case of no 


external magnetic field. As exhibited in Fig. 8 , the external 
magnetic flux is directed along the 2 : direction, normal to the 
junction face which resides in the xy plane so that the S/F 
and F/F interfaces are parallel to the y axis. We consider a 
rather general situation for the lengths and magnetization di¬ 
rections of the F strips. The lengths of the F strips can be 
unequal and are labeled by dpi and dp 2 (dpi 7 ^ dp 2 ). The 
strength of the exchange fields in both F regions are equal 
1 ^ 1 1 = 1 ^ 2 ! = 5 .OA 0 , while their orientations take arbitrary 
directions hi ^2 = (^T 2 : 21 2 )- Also, the junction width 

is equal to the width of the F stripes, i.e., Wpi = Wp 2 = 
Wp. 


The charge supercurrent given by Eq. (37) is comprised of 
even-frequency singlet terms, S, and the triplet odd-frequency 
components, T. To study exactly the behavior of each compo¬ 
nent of the Josephson charge current, we introduce the follow¬ 
ing decomposition scheme which is based on the discussions 
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FIG. 10. Spatial maps of the critical current density and pair potential as well as the critical charge current components for a junction 
configuration shown in Fig. 8. The same investigation as Fig. 9(b) is done here except we set dpi = dF 2 = l.O^s. The junction width and 
external magnetic flux are Wp = 2.05s and '1’ = S'l’o^ respectively. The vertical lines in the panels for the critical charge current components 
show the spatial separation of the two magnetic regions labeled Fi and F 2 . The quantities jf (a:, y) show the current components for flow 
along the x direction whereas Jy{x, y) represents the components for flow along the y direction. Here P denotes a specific component of the 
decomposed charge supercurrent namely, sO, sx, sy, and sz. 


and notations in the previous sections: 


terms involving S 


TSO 
^ x,y^ 

(39a) 

terms involving T^, 


TSX 

^x^y 5 

(39b) 

terms involving 


jsy 
^x^y 5 

(39c) 

terms involving 


JSZ 

^x,y’ 

(39d) 


In order to explicitly study the influence of each component of 
the decomposed charge supercurrent, we fix the magnetization 
of the Fi wire to be in the y-direction: hi — (0, , 0). Simi¬ 

larly, the magnetization of F 2 is orthogonal to that of Fi, and 
oriented in the z direction: /12 = (0, 0,h^). This orthogonal 
magnetization configuration is an inhomogeneous magnetic 
state that results in effective generation of equal-spin triplet 
correlations^. 

We show the spatial behavior of the total charge super¬ 
current density, pair potential, and the supercurrent compo¬ 
nents for a phase difference of |(^| = 7r/2 in Fig. 9 for the 
SIF/F/S junction depicted in Fig. 8. The geometric dimen¬ 
sions correspond to dpi = 0 . 25 s, dF 2 = l-85s, the junc¬ 
tion width Wp is set equal to 2.05s- The selected values of 
dp — dpi + dp 2 and Wp are useful towards understand¬ 
ing and analyzing the spin-triplet correlations in the different 
SI Ft FIS Josephson junctions considered here. For the panels 
found in (a), the external magnetic flux is absent, while for 
those in (b), a magnetic flux $ = 3$o is applied to the sys¬ 
tem. In panels (a), the charge current density is conserved as 
exhibited by its spatial uniformity throughout the ferromag¬ 
netic wire regions. There is also no y variations due the lay¬ 
ered system exhibiting translational invariance in that direc¬ 
tion. Without an external magnetic field, the charge supercur¬ 
rent (and pair potential) thus have no components along the 
y direction. The pair potential is shown to be an asymmetric 
function of position along the junction, which is simply due to 
the unequal length of each F wire and their differing magne¬ 


tizations. The charge supercurrent density components, 

and are also shown. These quantities can only vary 
spatially in the x-direction too since in the absence of an ex¬ 
ternal field, the total supercurrent is generated from the phase 
difference (j) between the S terminals, which vary only along 
the junction length {x direction). The odd-frequency triplet 
components of the supercurrent are localized in the Fi or F 2 
regions depending on the specific magnetization orientation in 
each region (see also the discussion in Sec. III). It is evident 
that disappears in the right F while is zero inside 
Fi segment where hi = { 0 ,h\, 0 ),h 2 = (0,0, h|), respec¬ 
tively (hi —)• J®^, h 2 —>■ T®^). Thus J®*' is generated in Fi 
and at the FIF interface is converted into J®^ in F 2 , or vice 
versa. The two components are generated in one F and do 
not penetrate into the other F with orthogonal magnetization. 
Turning now to Fig. 9(b), a finite magnetic flux of fl) = 3<l>o 
results now in a nonuniform supercurrent response that varies 
in both the x and y directions. Examining now the singlet cor¬ 
relations, the pair potential asymmetry is again present due to 
the spatially asymmetric magnetization regions, while the ex¬ 
ternal magnetic field induces vortices with normal state cores 
along the junction width. The singlet and triplet components 
of the supercurrent are plotted as well. The external magnetic 
field makes the problem effectively two dimensional, as it in¬ 
duces nonzero current densities in the y direction. The am¬ 
plitude of the equal-spin triplet component in the a;-direction, 
JJ^(a;, y), is somewhat smaller than the other components due 
in part to the small Fi width, dpi = 0.25s. However, compar¬ 
ing J^y{x,y) in panels (a) and (b), we see that the presence 
of a magnetic field can result in this triplet component now 
existing in F 2 , and it can also be more prominent in Fi. 

Turning now to the y components of the currents, we 
see that the magnitudes have increased in most instances 
by almost an order of magnitude or more. Although 
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FIG. 11. Spatial maps of the maximum total charge current, pair potential, and the even- {J^%{x,y)) and odd-frequency 
(Jx%{x,y), Jx%{x,y), Jx%{x,y)) components of the charge current. The configuration of the S/FIF/S Josephson junction considered 
here is shown in Fig. 12. The F/F interface is now perpendicular to the S/F interfaces. The external magnetic flux is zero, T* = 0, and the 
junction length dp is equal to 2.0^s. The width of each ferromagnetic strip are unequal, i.e., Wfi = 0.2^3 and Wf2 = l.S^g. The current 
components are plotted against the lateral y coordinate, i.e., along the junction width. Four various locations along the junction length are 
considered: x = 0.25^s, 0.50^s, 0.75^s, and l.O^s. Vertical lines in the panels for the current components separate the two ferromagnetic 
regions along the y direction (see Fig. 12). 


Jy°{x,y), Jyy{x,y), and Jy^{x,y) vary little at the differ¬ 
ent y locations, the overall spatial behavior is different for 
the odd and even frequency triplet components: The equal 
spin Jyy{x,y) penetrates considerably more into the ferro¬ 
magnetic regions compared with its x-component counterpart 
and compared with the spin-0 triplet Jy^{x, y). These results 
are consistent with the generation of a Meissner supercurrent 
in S/F/FIN structures. 

To study the effects of varying the system geometry and 
its corresponding effects on the triplet supercurrents, we con¬ 
sider more symmetric F regions with equal lengths dpi = 
dF 2 = l O^s. All other parameters are identical to the pre¬ 
vious S/F/FIS junction in Fig. 9(b). The first two panels on 
the left of Fig. 10 illustrate 2D spatial maps of the total su¬ 
percurrent density and pair potential. The remaining panels 
contain the spatial behavior of the singlet and triplet contri¬ 
butions to the total current. The effect of equal layer widths 
is reflected in the regular array of vortex patterns and circu¬ 
lating currents in the junction. The pair potential at x = 0 
exhibits three zeros when $ = 3$o. This is in accordance 
with the analytic expression in Eq. (30) and Eq. (38) which 
demonstrates the actual behavior of singlet correlations, con¬ 
stituting the spatial profile of proximity pair potential. To re¬ 
veal the even and odd frequency contributions to the contour 
plots, we examine in the remaining panels the supercurrent 
components as a function of x. The vertical lines separate the 
two ferromagnetic regions labeled by Fi and F 2 , as shown 
in the schematic of Eig. 8. We see that the formerly long- 
ranged Jyy{x,y) and J®*(x, j/), in Eig. 9 now vanish here 
when dpi = dp 2 = l-O^s. This demonstrates the crucial role 
that geometry can play in these type of junctions, in particu¬ 
lar the existence of odd-frequency correlations tend to favor 
configurations where the length of the ferromagnetic layers 
are unequal, for instance, dF 2 ^ dpi can further enhance 


the effect. This finding is also consistent with the results of 
a SIF!FIN configuration subjected to an external magnetic 
field,"^"^ where the component of the Meissner current is 
optimally anomalous if fiF 2 ^ dpi. The behavior of the sin¬ 
glet and triplet correlations discussed here are typically highly 
dependent on the magnetization of the double F layer, as well 
as the presence of an external magnetic field. In the next 
section, we therefore proceed to investigate another type of 
SIFIFIS Josephson junction with a different double F layer 
arrangement. 



FIG. 12. Schematic of the SIFIFIS Josephson junction where the 
interface of the double ferromagnet region is oriented perpendicular 
to the SIF interfaces. The junction is located in the xy plane and the 
SIF interfaces are along the y axis while the FIF interface is along 
the X axis. The width of the rectangular F strips Wfi, and Wf2 are 
generally unequal, i.e., Wfi 7^ Wf2- The lengths of the F strips, 
dp, however, are the same and correspond to the junction length. 
An external magnetic field H is applied normal to the junction plane 
along the 2 axis. 
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FIG. 13. Spatial maps depicting a vector plot of the total charge current density and density plot of the pair potential for the system shown in 
Fig. 12. The remaining line plots correspond to the charge supercurrent components: the singlet, y), and triplets, y), Jx%{x, y). 

Here the system parameters are identical to those used previously in Fig. 11, except the width of the ferromagnetic strips here are equal: 
Wfi = Wf2 = l.O^s. The junction length (If is equal to 2.0^s and there is no external magnetic flux, ih = 0. The vertical lines indicate the 
locations of the F/F interfaces in between the two ferromagnetic regions Fi and F 2 . 


IV.5. SIFIFIS heterostructures: F/F interface perpendicular 
to S/F Interfaces 

Here we consider a type of S/F/F/S Josephson junction 
where the previous ferromagnet double layer system is rotated 
by 7 r/ 2 , so that the F/F interface is normal to the outer S/F 
interfaces. This layout is shown in Fig. 12. As before, we 
take the length of the ferromagnetic strips to be equal, dpi — 
dF 2 = dp- We first consider ferromagnets where the width of 
the Fi and F 2 regions are unequal: Wpi ^ Wp 2 . When an 
applied magnetic field, FI, is present, it is directed along the z- 
axis, normal to the plane of the system. To be consistent with 
the previous subsection, we study two regimes of ferromag¬ 
netic widths namely, Wpi ^ Wp 2 and Wpi = Wp 2 , and 
investigate how the results vary for a finite external magnetic 
flux $. 

In Fig. 11, we show the results for Wpi — 0.2^s, and 
Wp 2 = 1 . 8 ^s, in the absence of an applied magnetic held 
($ = 0). First, the spatial map of the total maximum charge 
supercurrent shows that the current near the F/F junction has 
a non-zero y-component. This is in contrast to the uniform 
supercurrent shown in Fig. 9(a) for a “parallel” F/F segment. 
The induced y component is present over the entire junction 
width as exhibited by the distorted quasiparticle paths. The 
pair potential shows a symmetric behavior along the junction 
length (the x direction). However, similar to the previous 
case, the pair potential is an asymmetric function of y coor¬ 
dinate. This again arises from the unequal ferromagnetic wire 
widths (Wpi 7 ^ Wp 2 ) and their different magnetization ori¬ 
entations. The spatial behavior of the triplet and singlet con¬ 
tributions to the supercurrent is also shown as a function of 
position along the junction width (the y direction) at four x 
locations: x = 0.25^s, 0.5^s, 0.75^s, l.O^S. The amplitude 
of is around « 1/10 of and Jy^ but comparable with 
\Jyy -F Jy^\- These components of Jy, and hence Jy itself, 


vanish at y = ±0.5, corresponding to the vacuum boundary, 
and consistent with the boundary conditions given by Eq. (4). 
The vector plot of the charge supercurrent in Fig. 11 reveals 
no current flow along y in the middle of the junction, a; = 0 , 
throughout the junction width. The x component, is 
most negative at a: = 0 along the y axis while and 
are maximal. Therefore, the singlet-triplet conversion with 
Wp 2 Wpi is maximal at a: = 0. Another important as¬ 
pect of this S/F/F/S junction is seen in the behavior of 
and as a function of y: the two components are gener¬ 
ated in one F region and penetrate deeply into the adjacent F 
segment. By comparing these plots with those of a S/F/F/S 
configuration where F/F interface is parallel with the S/F in¬ 
terface subjected to an external magnetic held, presented in 
Fig. 9(b), one concludes that the spin-1 triplet components of 
charge supercurrent can penetrate the ferromagnetic regions 
when they flow parallel to the junction interfaces. 

In Fig. 13, we investigate the geometrical effects on the su¬ 
perconducting properties for a system with Fi and F 2 of equal 
width. We have Wpi = Wp 2 = l.O^s, and all previous pa¬ 
rameters remain intact. The current density map shows that 
the the induced y-component is equally distributed with re¬ 
spect to y = 0, the location of the F/F interface. The pair po¬ 
tential is a symmetric function of both the x and y coordinates 
about the origin, reflecting the symmetric geometric configu¬ 
ration. The panels representing the current components reveal 
that the y-component to the current density originates mainly 
from the singlet term, since it is evident that the odd- 

frequency components and are of opposite sign so 

that Jy^ ± J®^ « 0. Note that is an odd function of y 
for each fixed x-location. This reflects the fact that the junc¬ 
tion is symmetric along the y direction (see Fig. 11). From 
the current density spatial map the same functionality appears 
for the X coordinate. The integration dxJy^ gives the 
total charge supercurrent flowing in the y direction which is 
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FIG. 14. Spatial maps of total critical charge current density and pair potential and the decomposed components of charge supercurrent 
y), Jx^y{x, y), Jx^y{x, y), Jx%{x, y). The same parameter set as Fig. 11 are used except we now assume the system is subject to an 
external magnetic flux = S'Fo. The vertical lines in the current components’ panels separate the two ferromagnetic segments, labeled by Fi, 
and F 2 , with orthogonal magnetization orientations. The current components are plotted as a function of positions along the junction width 
i.e. y direction at four locations along the junction length i.e. x = 0.25^s, O.SO^s, 0.75^s, and l.O^s. 


therefore zero due since is antisymmetric in y. For the 
previous case shown in Fig. 11 , both triplet components 
and contributed to the induced y component in the total 
current density. Considering now the flow of current along 
the X direction, the top panels display the behavior of the odd- 
frequency triplet components with spin-1 and spin-0 projec¬ 
tions on the z quantization axes and respectively). 
Interestingly, the generation of each of the two triplet compo¬ 
nents tends to mirror each others behavior. Since the ferro¬ 
magnetic regions are symmetric, the singlet-triplet conversion 
of charge supercurrent components at the middle of junction 
{x = 0) across the junction is not as extensive compared to 
the previous case where WV 2 3> Wpx. It is clear that the total 
current passing through the junction must be conserved along 
the junction length (Itot{x) = const.). However, the singlet- 
triplet conversion in the charge supercurrent density is maxi¬ 
mal at a; = 0 and increases as the condition Wp 2 S> Wpi is 
fulfilled. 

In the previous two figures, the source of the driving current 
was the macroscopic phase differences between the S termi¬ 
nals in the Josephson structures. We now introduce an addi¬ 
tional source and associated triplet correlations by applying an 
external magnetic field (with corresponding flux $ = 3<I>o)^ 
normal to the junction plane along the z-axis (see schematic. 
Fig. 12). Therefore, Fig. 14 exhibits our results for a junction 
with the same parameters used in Fig. 11 but now, $ = 3<l>o. 
The magnetic field causes the quasiparticles to undergo cir¬ 
cular motion as shown in the vector plot for the supercurrent 
density. The pair potential also vanishes at particular loca¬ 
tions, which now form a nontrivial pattern. The charge su¬ 
percurrent components are clearly modified by the external 
magnetic field compared to $ = 0 (Fig. 11). We have found 
that the behavior of the supercurrent components along y are 
similar to the corresponding components along x when the 
FIF junction is parallel with the SIF interfaces (Fig. 9(b)). 
In both cases, the supercurrent flows normal to the FIF junc¬ 


tion. Likewise, the component as a function of y in Fig. 11 
is similar to the Jy component vs. x as seen in Fig. 9(b). Re¬ 
call that in Fig. 14, we have fl) = 3$0 j while in Fig. 9(b) 
$ = 0. Here we And also that the triplet components of the 
charge supercurrent generated in one F strip, flowing parallel 
with the FIF interface deeply penetrate the adjacent F com¬ 
pared to when the current flows normal to the FIF interface. 
It is also evident that the total net supercurrent in the y direc¬ 
tion is zero, i.e., dx{Jy^ + Jy^ + Jy^) = 0. 


Upon changing the width of the F layers toWpi = Wp 2 = 
l.O^S, Fig. 15 shows that the circular paths observed in the 
current density and zeroes in the spatial map of the pair po¬ 
tential have reverted back to a more symmetric configura¬ 
tion compared to Fig. 14. Here the magnetization orienta¬ 
tions in the F regions are orthogonal: hi = (0,and 
^2 = (0, 0,/i^). The panels containing the current compo¬ 
nents and as a function of y are now much more lo¬ 
calized in the F regions, in contrast to the previous case where 
Wp 2 S> Wpi (Fig. 11). Considering the current density along 
the y direction, the triplet components are also localized in the 
F regions. Similarities are observed when comparing their y 
dependence with the x dependence of the triplets and J®^ 
in Fig. 10. This comprehensive investigation into the different 
SI FI FIS structures has shown that an applied magnetic held 
can result in the appearance of odd-frequency triplet compo¬ 
nents to the charge transport when the supercurrent is paral¬ 
lel to the FIF interface. These odd-frequency correlations 
can exhibit extensive penetration into ferromagnetic regions 
with orthogonal magnetizations. These findings are in stark 
contrast to those cases where the supercurrent flows normal 
to uniform FIF double layers with orthogonal magnetization 
orientations (see also the discussions of one-dimensional sys¬ 
tems in Sec. III). 
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FIG. 15. Spatial maps of the total charge supercurrent and pair potential in addition to the current components 

Jx%{x,y), Jx^y{x,y), Jx^y{x,y), J^^y{x,y), at maximum charge current flowing across the junction. The same investigations as Fig. 13 
are presented here. The junction length is fixed at = 2.0^s while the width of each ferromagnetic strip, labeled by Tj and F 2 , are 
Wfi = Wf2 = l.O^s. Unlike Fig. 13, the external magnetic flux is 'F = S'Fq. The two ferromagnetic regions with orthogonal magnetization 
orientations are separated by vertical lines and labeled Fi and F 2 . 


IV.6. Spin valve structure probing of equal-spin triplet 
supercurrent 


Having analyzed in detail the singlet-triplet contributions 
to the supercurrent in various situations, we now proceed to 
demonstrate an experimentally accessible structure that can 
directly detect our predictions involving parallel transport in 
F/F bilayers. We convert one of the outer S terminals to 
a finite sized normal layer, so that in effect we consider a 
simple S/F/F/N spin valve (Fig. 16). This structure gener¬ 
ates pure odd-frequency spin-1 triplet correlations, and can be 
a more experimentally accessible system for generating and 
controlling triplet correlations in supercurrent transport.^ The 
basic structure is made of two uniform ferromagnetic layers 
with thicknesses dpi, and dF 2 , and a relatively thick nor¬ 
mal layer dN, which is connected to a superconducting ter¬ 
minal. The N layer assists in probing the equal-spin compo¬ 
nent of the supercurrent efficiently since the exchange field 
there is zero, and this component of supercurrent decays very 
slowly in the normal metal. The magnetization of the left 
F is determined by a rotation angle /3i with respect to the 
z axis {hi — /io(cos7r/4sin/3i,sin7r/4sin/3i,cos/3i)) while 
the magnetization of right F layer is assumed to be fixed along 
the z axis. The non-superconducting part of the spin-valve 
is subject to an external magnetic field H oriented along the 
z axis, Hz- The external magnetic field leads to a diamag¬ 
netic supercurrent J, which depends on x and flows along the 
y direction [denoted by Jy{x)'\ parallel to the interfaces. In 
the presence of an external magnetic held, we make the usual 
substitution 5 —> V + eA/c, where A is the vector potential, 
related to the external magnetic held (i?) via V x A = H. In 
the linear response regime^^^®®, the supercurrent density can 


be expressed by: 

Jy{x) = —Jo8ieAy{x) x 

(ietanhe/3 {S(e, a:)S*(—e, x) — S(—e, x)S*(e, x)+ 

T 3 ;(-e, a;)T*(e, x) - T 3 ;(£, a;)T* (-e, x)+ 

Ty(-£:, x)T*{e, x) - Ty(e, a:)T* (-e, a:)-F 
Tz{-e,x)Tl[e,x) - Tz(e,x)Tl{-£,x)} . (40) 

Strictly speaking, to determine Ay{x), it is necessary to 
solve Maxwell’s equation incorporating the Coulomb gauge 
in conjunction with the expression for Jy{x) and appropri¬ 
ate boundary conditions for We have found however 




h={oaha 


FIG. 16. (Color online) Schematic of the proposed S/F/F/N junc¬ 
tion. The magnetization orientation of the middle F layer with 
width dF 2 is fixed along the z direction. The left F layer, with 
width dFi, has a misaligned magnetization direction correspond¬ 
ing to hi — fio(cos7r/4, sin7r/4sin/3i, cos/3i). The normal metal 
layer has width dpf. An external magnetic field is applied to the sys¬ 
tem in the z direction parallel to the interfaces: H = (0, 0, Flz). Due 
to the external magnetic field, a diamagnetic supercurrent flows in the 
y direction J{x) — (0, Jy{x), 0) and varies along x. The junction 
interfaces lie in the yz plane. 
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FIG. 17. (Color online) Top row. Spatial behavior of the total supercurrent and its components as a function of position inside the S/F/F/N 
junction for various magnetization orientations in Fi: hi = /io(cos7r/4sin/3i, sin7r/4sin/?i, cos/3i) for/3i = 0 , 0 .lvr, 0 . 27 r, O.Stt, 0.47r, 
and O.Stt. The thickness of the left F layer is fixed at dpi = while the thickness of the right F layer is equal to d_F 2 = l.SS^g. 

Bottom row. The left and right thicknesses in the F layers are now equal to dpi = dp 2 = is- The sum of dpi and dp 2 is identical for both 
cases investigated in the top and bottom rows. The normal layer thickness is also fixed at 2.5^s in both cases. The logarithmic scales on the 
x-axes permits easier visualization of the long-ranged spatial behavior of the odd-triplet components that have spin projections normal to the 
z quantization axis. 


that through our extensive numerical investigations, Ay{x) 
is typically a linear function of x and weakly varies with 
the magnetization alignment, /3i. Therefore, it is the en¬ 
ergy and spatial dependence of the Green’s function com¬ 
ponents in Eq. (40) which governs the supercurrent density. 
To study the supercurrent behavior, and similar to what was 
done above for Josephson structures, we separate out the 
even- and odd-frequency components of the total supercur¬ 
rent the same as Sec. III.l. We achieve this by first defin¬ 
ing Itot{x)=Iso{x)+Isx{x)+Isy{x)+Isz{x), where Iso{x), 
Isx{x), Isy{x), and are given by, 

Iso{x) = 

de {S(e)S*(-e) - S(-e)S*(e)} tanhe/3, (41a) 

d^Sxi^X^ — 

de {Ta;(—e)T*(e) — Ta;(£)T* (—e)} tanhe^,(41b) 

Isy(.x) = 

de {Ty(-e)T*(e) - Tj^(e)T*(-£)} tanhe/3,(41c) 
Iszix) = 

de {T^(-e)T;(e) - T^(e)T;(-e)} tanhe/3. (41d) 

Also the boundary conditions at the right far end interface 
of the normal metal should be modified accordingly (see 

Fig. 12 ): 


dx{l-'^xi~£) + iTy^—e)) — 0 , 

(42a) 

9,(TS(-e)+T,(-e)) =0, 

(42b) 

9,(TT:(e)-*T;(e))=0, 

(42c) 

5,(TS*(e)+T:(e))=0. 

(42d) 


Having now established the method in which to determine 
the spatial behavior of the supercurrent and its decomposi¬ 






tions, we present in Fig. 17 our findings for the proposed 
nanovalve subject to an external magnetic field. Several val¬ 
ues of the left magnetization orientation /3i are investigated 
(the right F magnetization is fixed along the z axis), corre¬ 
sponding to Pi/n = 0,0.1, 0.2,0.3, 0.4,0.5. In the top row, 
the thickness of the left F layer is much smaller than the 
thickness of the adjacent F layer: dFi{= 0.15^s) dF 2 {= 

1.85^s). The thickness of the normal metal is of moderate 
size, set equal to dM = 2.5^s. The total non-superconducting 
region thus has a thickness of 4.5^s. To demonstrate the long- 
ranged nature of the spin -1 triplet components Isx{x) and 
Isy {x) in this simple structure with only two magnetic regions 
with uniform and non-collinear magnetization orientations, 
we introduce a three dimensional magnetization in the left F 
region: hi — /io(cos 7 r/ 4 sin/ 3 i, sin 7 r/ 4 sin, 5 i, cos/3i). To 
clearly show the spatially long-ranged components, we use a 
logarithmic scale for each x-axis. It is seen that the triplet 
component with zero spin projection on the z-axis, Isz{x), 
is short-ranged and becomes suppressed at positions approx¬ 
imately greater than x ^ 0.5^s. This is in drastic contrast 
to the Isx{x), and Isy{x) components, which propagate over 
the entire non-superconducting region. The singlet compo¬ 
nent, Iso{x), vanishes around x ^ which is about half the 
thickness of the magnetic layers (d^i + dF 2 = 2 ^ 5 ). Thus, 
the long-ranged components can strongly influence the total 
supercurrent, Itot{x), and be the primary source of nonzero 
supercurrent in the normal metal. These findings are there¬ 
fore suggestive of an experimentally accessible S/F/F/N 
spin valve with long-ranged triplet correlations controllable 
by magnetization rotation in one of the F layers. 

It turns out that the the pure spin-1 odd-frequency spin- 
valve effect is optimal for asymmetric F layer thicknesses. 
To illustrate this, we consider a symmetric configuration with 
dFi — dF 2 = is, which ensures the same total thickness 
of the ferromagnet regions (the normal metal thickness is un¬ 
changed at dN = 2.55s). The results of this configuration 
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are represented in the bottom row of Fig. 17. The total cur¬ 
rent and its components are plotted as a function of position 
inside the non-superconducting regions for the same magneti¬ 
zation orientations in the asymmetric case (top row). We find 
that in contrast to the panels in the top row, the triplet compo¬ 
nents Isx{x), Isy{x), and Iszix) are now all short-ranged in 
the sense that they vanish near x ~ 0 . 5 ^ 5 . In particular, the 
spin -1 components, Isxix), and Isy{x) do not propagate in 
F 2 , where the magnetization is orthogonal to the correlations’ 
spin orientations. The singlet contribution, Iso{x), for both 
the dp 2 ^ dpi and dp 2 = dpi cases shows similar behavior, 
where it vanishes at around x ~ ^s- This can be understood 
by noting that the spin-splitting effects of ferromagnetism de¬ 
stroy opposite-spin paring correlations. The total current also 
vanishes at x ~ ^5 for all magnetization orientations, /3i, in 
contrast to what was observed in the dp 2 ^ dpi structure 
above. 

We are now able to compare the behavior of various su¬ 
percurrent components in the S/FI FIN and the SI FI FIS 
Josephson structures presented in Sec. III. As seen, our in¬ 
vestigations for a wide range of parameters, including mis¬ 
aligned SIFIFIS structures with dp 2 ^ dpi, illustrate local¬ 
ized spin -1 components, Isx{x), Isy(x), whereby their spin 
orientation follows the local magnetization orientation. These 
components of supercurrent are unable to propagate through¬ 
out the bilayer FIF junctions where there is a relative orthog¬ 
onal magnetization. This comparison immediately reveals the 
efficiency and advantages of a SIFIFIN nanovalve for de¬ 
tecting proximity spin -1 triplet superconducting correlations 
experimentally. The pure spin-1 triplet diamagnetic response 
discussed here for the SIFIFIN spin valves might be mea¬ 
sured using a local probe such as the density of states or cur¬ 
rent density. As is well understood, the triplet correlations 
can alter the induced minigap in the local density of states in 
the nonsuperconducting regions. For example, the traditional 
s-wave minigap profile can become peaked due to the emer¬ 
gence of triplet correlations which cause a resonance near the 
Fermi surface. Also, as seen, our findings demonstrate that 
only spin -1 triplet components of the diamagnetic current sur¬ 
vive in the N region if dp 2 dpi, for sufficient magne¬ 
tization misalignment. Therefore, these results suggest that 
currents measured experimentally within the N wire for such 
systems should be comprised of purely equal-spin triplet cor¬ 
relations that are odd-frequency in character. 

Therefore, the spin-parameterized technique reveals the 
spatial profile of the odd and even frequency components of 
the total supercurrent and provides additional insight into their 
behaviors. The general three-dimensional approach which we 
outlined can then be employed for two-dimensional finite-size 
magnetic/superconducting proximity hybrids with arbitrary 
magnetization patterns. In this paper we categorized SIF hy¬ 
brids with layered F regions into two classes: configurations 
hosting supercurrent transport 1 ) parallel and 2 ) perpendicular 
to the FIF interface. The one-dimensional structures studied 
in Sec. Ill belong to class 2), while the two-dimensional con¬ 
figurations discussed in this section are examples of class 1 ). 
For the fairly wide range of configurations considered in this 
paper, we demonstrated that singlet supercurrent flow “par¬ 


allel” to uniformly magnetized F/F interfaces can generate 
substantial long-range triplet supercurrents regardless of the 
specific configuration and geometry^’'®’^*. This is of course in 
stark contrast to their counterparts where the singlet supercur¬ 
rent flows “perpendicular” to the FIF cross sections and fully 
discussed in Sec. III. 


V. CONCLUSION 

In conclusion, we have employed a quasiclassical method 
and spin-parameterization technique in both the Usadel equa¬ 
tions and associated boundary conditions. This approach pro¬ 
vides a suitable computational and analytical framework to 
pinpoint the spatial behaviors of the decomposed odd and 
even frequency supercurrent components in layered ferro¬ 
magnetic junctions. We have studied the transport characteris¬ 
tics of the spin decomposed supercurrent components for two 
generic scenarios and several experimentally relevant struc¬ 
tures: systems with 1 ) parallel or 2 ) perpendicular charge su¬ 
percurrent flow relative to the magnetic interfaces. Two types 
of finite-size two-dimensional magnetic SIFIFIS Josephson 
junctions subject to an external magnetic field, and support¬ 
ing supercurrents flow parallel to the FIF interface, are con¬ 
sidered. In one type, we assume the FIF interface is parallel 
to the SIF interfaces, while in the other type, the FIF in¬ 
terface is perpendicular to the SIF interfaces. Our studies 
revealed that when a supercurrent flows parallel to the FIF 
junction, the long-ranged spin-triplet supercurrent will be ef¬ 
fectively generated and it will propagate deeply into the fer¬ 
romagnetic regions independent of junction geometries. This 
phenomenon is more pronounced when the thickness of fer¬ 
romagnetic strips are unequal (dp 2 dpi or Wp 2 S> Wpi), 
depending on the junction type under consideration. This 
effect disappears when the supercurrent direction is perpen¬ 
dicular to the FIF interface. To gain more insight and to 
have explicit comparisons, we also studied the various aspects 
of the spin-decomposed supercurrent components in one¬ 
dimensional S/F/F/S, S/F/F/F/S, and SjHojFiHojS 
hybrids using the spin-parametrization scheme. We studied 
the spatial maps of each supercurrent component and found 
that the total supercurrent is dominated by the long-range 
triplet component, which also governs the behavior of total 
supercurrent. Our results demonstrated that structures with 
inhomogeneous magnetization patterns generate long-ranged 
supercurrents stronger than their uniform counterparts and 
that the component of supercurrent corresponding to the ro¬ 
tating magnetization component is long-ranged in such struc¬ 
tures. Finally, we proposed an SIFIFIN spin-valve, which 
represents an experimentally accessible platform to probe the 
predicted phenomena in the parallel supercurrent transport 
scenario. We have studied the spatial maps of the singlet 
and triplet supercurrent components of this spin-valve when 
subjected to an external magnetic field. Our findings demon¬ 
strated that for misaligned magnetizations and different thick¬ 
nesses of the F layers, dp 2 dpi, a long-ranged odd-triplet 
component of the supercurrent response arises and populates 
the adjacent normal metal N. 


24 


ACKNOWLEDGMENTS 

We would like to thank G. Sewell for his valuable instruc¬ 
tions in the numerical parts of this work. We also appreciate 
N.O. Birge and F.S. Bergeret for useful discussions and com¬ 
ments. K.H. is supported in part by ONR and by a grant of 
supercomputer resources provided by the DOD HPCMR 


* phymalidoust@gmail.com 
^ klaus.halterman@navy.mil 

* I. Zutic, J. Fabian, and S.D. Sarma, Rev. Mod. Phys.76, 323 
(2004). 

^ A.A. Golubov, M.Yu. Kupriyanov, E. Ilichev, Rev. Mod. Phys. 
76,411 (2004). 

^ A. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 

F. S. Bergeret, A.F. Volkov, K.B. Efetov, Rev. Mod. Phys. 77, 
1321 (2005). 

® K.B. Efetov, I.A. Garifullin, A.E. Volkov,K. Westerholt, Mag¬ 
netic Heterostructures,Advances and Perspectives in Spinstruc- 
tures and Spintransport. ed. by H. Zabel, S.D. Bader, Series. 
Springer Tracts in Mod. Phys. vol 227 (Springer, New York, 
2007), P. 252. 

® M. Eschrig, Phys. Today 64, 43 (2011). 

’ M. Alidoust and K. Halterman, Appl. Phys. Lett. 105, 202601 
(2014). 

* M. Alidoust and K. Halterman, Phys. Rev. B 89, 195111 (2014). 
’ K. Halterman, O. T. Vails, and M. Alidoust, Phys. Rev. Lett. Ill, 

046602 1694 (2013). 

A.M. Bobkov and I.V. Bobkova, Phys. Rev. B 84, 054533 (2011); 
Phys. Rev. Lett. 108, 197002 (2012). 

" A.G. Malshukov, and A. Brataas, Phys. Rev. B86, 094517 (2012). 
L.B. Ioffe, V.B. Geshkenbein, M.V Feigelman, A.L. Fauchere, 
and G. Blatter, Nature (London) 398, 679 (1999). 

Y. S. Barash, I.V. Bobkova, and T. Kopp, Phys. Rev. B 66, 
140503(R) (2002). 

A. Cottet, T. Kontos, W. Belzig, C. Schonenberger, C. Bruder, 
Europhys. Lett. 74, 320 (2006). 

B. Crouzy, S. Tollis, and D.A. Ivanov, Phys. Rev. B 76, 134502 
(2007). 

Ya.V. Fominov, A.F. Volkov, and K.B. Efetov, Phys. Rev. B 75, 
104509 (2007). 

Z. Pajovic, M. Bozovic, Z. Radovic, J. Cayssol, and A. Buzdin, 
Phys. Rev. B 74, 184509 (2006). 

H. Sellier, C. Baraduc, F. Lefloch, and R. Calemczuk, Phys. Rev. 
Lett. 92, 257005 (2004). 

N.G. Pugach, M.Yu. Kupriyanov, A.V. Vedyayev, C. Lacroix, E. 
Goldobin, D. Koelle, R. Kleiner, and A.S. Sidorenko Phys. Rev. 
B 80, 134516(2009). 

L. -J. Jin, Y. Wang, L. Wen, G.-Q. Zha, and S.-P. Zhou Physics 
Letters A 376, 2435 (2012). 

M. Eschrig, J. Kopu, J.C. Cuevas, and G. Schon, Phys. Rev. Lett. 
90, 137003 (2003). 

M.N. Baibich, J.M. Broto, A. Fert, F. Nguyen Van Dau, F. 
Petroff, P. Eitenne, G. Creuzet, A. Friederich, and J. Chazelas, 
Phys. Rev. Lett. 61, 2472 (1988). 

G. Binasch, P. Grunberg, F. Saurenbach, and W. Zinn, Phys. Rev. 
B 39, 4828 (1989). 

M. Alidoust, G. Sewell, and J. Linder Phys. Rev. B 85, 144520 

( 2012 ). 


M. Alidoust, J. Linder, G. Rashedi, T. Yokoyama, and A. Sodbo, 
Phys. Rev. B81, 014512 (2010). 

V.I. Zdravkov, J. Kehrle, G. Obermeier, D. Lenk, H.-A. Krug von 
Nidda, C. Muller, M.Yu. Kupriyanov, A.S. Sidorenko, S. Horn, 

R. Tidecks, and L.R. Tagirov, Phys. Rev. B 87, 144507 (2013). 

S. V. Bakurskiy, N.V. Klenov, I.I. Soloviev, M.Yu. Kupriyanov, 
and A.A. Golubov, Phys. Rev. B 88, 144519 (2013). 

Y. Makhlin, G. Schoen, A. Shnirman, Rev. Mod. Phys. 73, 357 

( 2001 ). 

V.V Ryazanov, V.A. Oboznov, A.Yu. Rusanov, A.V. Vereten¬ 
nikov, A.A. Golubov, and J. Aarts, Phys. Rev. Lett. 86, 2427 
( 2001 ). 

L. N. Bulaevskii, M.J. Graf, C.D. Batista, and V.G. Kogan, Phys. 
Rev. B 83, 144526 (2011). 

A.I. Buzdin, L.N. Bulaevskii, S.V. Panyukov, JETP Lett. 35, 178 
(1982). 

K. Halterman, and O.T. Vails, Phys. Rev. B 69, 014517 (2004). 
F.S. Bergeret, A.F. Volkov, K.B. Efetov, Phys. Rev. Lett. 86, 4096 
( 2001 ). 

Y. Asano, Y. Tanaka, T. Yokoyama, and S. Kashi way a, Phys. Rev. 
B 74, 064507 (2006). 

M. Eschrig, T. Lofwander, Nat. Phys. 4, 138 (2008). 

T. Lofwander, T. Champel, J. Durst, M. Eschrig, Phys. Rev. Lett. 
95, 187003 (2005). 

T. Kontos, M. Aprili, J. Lesueur, X. Grison, Phys. Rev. Lett. 86, 
304 (2001). 

I. Sosnin, H. Cho, VT. Petrashov, and A.F. Volkov, Phys. Rev. 
Lett. 96, 157002 (2006). 

K. Halterman, PH. Barsic, and O.T. Vails Phys. Rev. Lett. 99, 
127002 (2007). 

T.S. Khaire, M.A. Khasawneh, W.P. Pratt Jr., N.O. Birge, Phys. 
Rev. Lett. 104, 137002 (2010). 

R. S. Keizer, S.T.B. Goennenwein, T.M. Klapwijk, G. Miao, G. 
Xiao and A. Gupta, Nature 439, 825 (2006). 

M. Houzet and A.I. Buzdin, Phys. Rev. B 76, 060504(R) (2007). 
M. Alidoust, and J. Linder, Phys. Rev. B82, 224504 (2010). 

** M. Alidoust, K. Halterman, J. Linder, Phys. Rev. B 89, 054508 
(2014). 

M. Houzet, A.I. Buzdin, Phys. Rev. B74, 214507 (2006). 

M. Knezevic, L. Trifunovic, Z. Radovic, Phys. Rev. B 85, 094517 

( 2012 ). 

L. Trifunovic, Phys. Rev. Lett. 107, 047001 (2011). 

T.Y Karminskaya, M.Y Kupriyanov, and A.A. Golubov, JETP 
Lett. 87, 570 (2008); A.I. Buzdin, A.S. Melnikov, and N.G. Pu¬ 
gach, Phys. Rev. B 83, 144515 (2011). 

S. Hikino and S. Yunoki, Phys. Rev. Lett. 110, 237003 (2013). 

J. W.A. Robinson, J.D.S. Witt, M.G. Blamire, Science 329, 5987 

( 2010 ). 

C. Wu, O.T. Vails, and K. Halterman Phys. Rev. Lett.108, 117005 

( 2012 ). 

K. Halterman and O. T. Vails Phys. Rev. B 80, 104502 (2009). 



25 


V.A. Oboznov, V.V. Bolginov, A.K. Feofanov, V.V. Ryazanov, 
and A.I. Buzdin, Phys. Rev. Lett. 96, 197003 (2006). 

G.A. Ovsyannikov, A.E. Sheyerman, A.V. Shadrin, Yu.V. Kislin- 
skii, K.Y. Constantinian, and A. Kalabukhov JETP Letters 97, 
145 (2013) 

T. Yu. Karminskaya, A.A. Golubov, and M.Yu. Kupriyanov, Phys. 
Rev. B 84, 064531 (2011). 

C. Wu, O.T. Vails, and K. Halterman, Phys. Rev. B 86, 014523 

( 2012 ). 

A. S. Melnikov, A.V. Samokhvalov, S.M. Kuznetsova, and A.I. 
Buzdin, Phys. Rev. Lett. 109, 237006 (2012). 

S. Oh, D. Youm and M.R. Beasley, Appl. Phys. Lett. 71, 2376 
(1997). 

YV. Eominov el al, JETP Lett. 91, 308 (2010). 

“ PV Leksin, N.N. Garifyanov, LA. Garifullin, Ya.V Eominov, J. 
Schumann, Y. Krupskaya, V. Kataev, O.G. Schmidt, and B. Bch- 
ner, Phys. Rev. Lett. 109, 057005 (2012). 

PV Leksin, N.N. Garifyanov, LA. Garifullin, J. Schumann, V. 
Kataev, O.G. Schmidt, and B. Bchner, Phys. Rev. Lett. 106, 
067005 (2011). 

PV Leksin, N.N. Garifyanov, LA. Garifullin, J. Schumann, V. 
Kataev, O.G. Schmidt, and B. Bchner, Phys. Rev. B 85, 024502 
( 2012 ). 

B. Li, N. Roschewsky, B.A. Assaf, M. Eich, M. Epstein-Martin, 

D. Heiman, M. Mnzenberg, and J. S. Moodera, Phys. Rev. Lett. 
110, 097001 (2013). 

^ PV Leksin, N.N. Garifyanov, LA. Garifullin, J. Schumann, H. 
Vinzelberg, V. Kataev, R. Klingeler, O.G. Schmidt and B. Bch¬ 
ner, Appl. Phys. Lett. 97, 102505 (2010). 

A.A. Jara, C. Safranski, I.N. Krivorotov, C. Wu, A.N. Malmi- 
Kakkada, O.T. Vails, and K. Halterman, Phys. Rev. B89, 184502 
(2014). 

“ M.G. Elokstra, T.C.Cunningham, J. Kim, N. Satchell, G. Bur¬ 
nell, S.J. Bending, P.J. Curran, S.J. Langridge, C. Kinane, J.F.K. 
Cooper, N. Pugach, M. Eschrig, S.L. Lee, arXiv: 1404.2950. 

V. Braude, Yu.V. Nazarov, Phys. Rev. Lett. 98, 077003 (2007). 

E. S. Bergeret and J.C. Cuevas, J. Low Temp. Phys., 153 304 
(2008). 

® J.M. Rowell, Phys. Rev. Lett. 11, 200 (1963). 

™ B.D. Josephson, Phys. Lett. 1, 251 (1962). 

J. C. Cuevas and F.S. Bergeret, Phys. Rev. Lett., 99 217002 
(2007). 

M. Alidoust, G. Sewell, and J. Linder, Phys. Rev. Lett. 108, 
037001 (2012). 

M. Alidoust, and J. Linder, Phys. Rev. B 87, 060503(R) (2013). 

U. Ledermann, A.L. Pauche‘re, and G. Blatter, Phys. Rev. B 59, 
R9027 (1999). 

K. Usadel, Phys. Rev. Lett. 25, 507 (1970); A.I. Larkin and YN. 
Ovchinnikov, in Nonequilibrium Superconductivity, edited by D. 


Langenberg and A. Larkin (Elsevier, Amesterdam, 1986), P. 493. 
J.P. Morten, M.S. thesis, Norwegian University of Science and 
Technology, 2003. 

W. Belzig, F.K. Wilhelm, C. Bruder, G. Schn, A. D.Zaikin, Su- 
perlattices and Microstructures, Vol. 25, Iss. 56, Pg. 1251 (1999). 
A. V. Zaitsev, Zh. Eksp. Teor. Eiz. 86, 1742 (1984) (Sov. Phys. 
JETP 59, 1015 (1984); M. Y. Kuprianov et al, Sov. Phys. JETP 
67, 1163 (1988). 

L. Angers, F. Chiodi, G. Montambaux, M. Ferrier, S. Gueron, H. 
Bouchiat, and J.C. Cuevas, Phys. Rev. B 77, 165408 (2008). 

F. Chiodi, M. Ferrier, Gueron, J.C. Cuevas, G. Montambaux, F. 
Fortuna, A. Kasumov, H. Bouchiat, Phys. Rev. B 86, 064510 
( 2012 ). 

J. R. Clem, Phys. Rev. B 81, 144515 (2010). 

A.F. Volkov and K.B. Efetov, Phys. Rev. B 81, 144522 (2010). 

T. Champel, M. Eschrig, Phys. Rev. B 72, 054523 (2005). 

M. Eschrig, T. Lofwander, T. Champel, J.C. Cuevas, J. Kopu, G. 
Schon, J. Low Temp. Phys. 147, 457, (2007). 

G. B. Halsz, J.W.A. Robinson, J.P. Annett, and M.G. Blamire, 
Phys. Rev. B 79, 224505 (2009) 

D. Fritsch and J.R Annett, New J. Phys. 16, 055005 (2014). 

H. Suderow, 1. Guillamon, J.G. Rodrigo, S. Vieira, 
arXiv:1403.5514. 

** A.A. Abrikosov, Fundamentals of the Theory of Metals (North- 
Holland, Amsterdam, 1988). 

T. Yamashita, L. Rindereer, J. Low Temp. Phys. 21, 153 (1975). 
A.A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 

M. Tinkham, Introduction to Superconductivity (McGraw-Hill, 
New York, 1996). 

T. Lofwander, T. Champel, M. Eschrig, Phys. Rev. B 75, 014512 
(2007). 

Higashitani and K. Nagai, Physica B 194, 1385 (1994). 

W. Belzig, C. Bruder, and G. Schon, Phys. Rev. B 53, 5727 
(1996). 

P. Visani et al., Phys. Rev. Lett. 65, 1514 (1990). 

F. Bemd, Muller-Allinger and A. Celia Mota, Phys. Rev. Lett. 84, 
3161 (2000). 

F. Bernd Muller-Allinger and A. Celia Mota, Phys. Rev. B 62, 
6120 (2000). 

C. Bruder and Y. Imry, Phys. Rev. Lett. 80, 5782 (1998). 

T. Yokoyama, Y. Tanaka, and N. Nagaosa, Phys. Rev. Lett. 106, 
246601 (2011). 

K. Halterman, O. T. Vails, and M. Alidoust, Phys. Rev. Lett. Ill, 
046602 1694 (2013). 



